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Um  das  vibroakustische  Übertragungsverhalten  von 
mechanischen  Systemen  rechnerisch  zu  bestimmen, 
sind  geeignete  Modelle  des  Realsystems  erforderlich. 
Nur  für  sehr  einfache  Systeme  kann  das  Systemver¬ 
halten  unmittelbar  mathematisch-analytisch  formu¬ 
liert  werden.  Für  komplexe  Systeme  sind  entspre¬ 
chend  vereinfachte  Modelle  erforderlich,  deren  An¬ 
wendungsbereich  aufgrund  der  getroffenen  Vereinfa¬ 
chungen  naturgemäß  begrenzt  ist.  Als  nichtkontinu¬ 
ierliche  Methode  für  niederfrequente  Übertragung  hat 
sich  die  Mehrkörpersimulation  (MKS)  etabliert.  Ty¬ 
pisch  sind  für  diese  Methode  die  Annahme  starrer 
Körper  sowie  lineare  und  nichtlineare  Koppelelemen¬ 
te.  Bei  räumlich  ausgedehnten  Strukturen  und  höhe¬ 
ren  Frequenzen  ist  die  Annahme  starrer  Körper  oft 
nicht  mehr  zulässig.  Hier  verwendet  man  die  Finite- 
Elemente-Methode  (FEM).  Dabei  wird  der  Körper  in 
viele  kleine  Elemente  unterteilt  (Quader,  Tetraeder, 
usw.).  Die  Stetigkeitsbedingungen  an  den  Begren¬ 
zung  sflächen  der  einzelnen  Flächen-  oder  Volumen¬ 
elemente  führen  letztlich  auf  ein  lineares  Gleichungs¬ 
system.  Die  Lösungen  des  homogenen  Gleichungs¬ 
systems  ergeben  die  Eigenschwingformen  und  Eigen¬ 
frequenzen  der  Struktur.  Im  höherfrequenten  Bereich 
nimmt  die  Anzahl  der  Schwingungsmoden  stark  zu, 
eine  exakte  Betrachtung  wird  dementsprechend  im¬ 
mer  aufwändiger.  Durch  statistische  Mittelung  über 
viele  Eigenschwingformen  kann  das  Problem  wieder 
vereinfacht  und  die  Rechenzeit  begrenzt  werden. 
Hier  setzt  die  statistische  Energieanalyse  (SEA)  an. 
Die  damit  verbundene  Näherung  ist  jedoch  nur  exakt 
genug,  wenn  die  Zahl  der  Resonanzen  entsprechend 
groß  und  damit  die  Frequenz  oberhalb  von  ca.  500  Hz 


Bild  15-1:  Gängige  Analyseverfahren  und  deren  typi¬ 
sche  Anwendungsbereiche 


liegt.  Während  sich  die  typischen  Frequenzbereiche 
von  MKS-  und  FEM-Berechnungen  weit  überlappen, 
setzt  die  SEA  erst  oberhalb  des  praktisch  nutzbaren 
Frequenzbereiches  der  FEM  von  maximal  ca.  200  Hz 
ein.  Man  spricht  daher  vom  „Mid-Frequency-Gap“ 
zwischen  200  Hz  und  500  Hz.  Dieses  durch  geeignete 
Methoden  zu  schließen  ist  derzeit  noch  Gegenstand 
der  Forschung. 

15.1  Mehrkörpersimulation  (MKS) 

Mehrkörpersimulation  findet  bei  der  Berechnung  des 
dynamischen  Systemverhaltens  in  der  Fahrzeugtech¬ 
nik  vielfältige  Anwendung.  Einsatzgebiete  sind  dabei 
sowohl  die  Fahrdynamik-Analyse  (Handling)  als 
auch  Komfortuntersuchungen  (Ride,  NVH).  Im  An¬ 
triebsbereich  werden  z.B.  der  Triebstrang  und  der 
Ventiltrieb  mit  Hilfe  von  MKS- Analysen  ausgelegt. 


Bild  15-2:  MKS  Modell  von  Triebstrang  und  Fahr¬ 
werk  mit  Hydropulsanregung 

Ein  dynamisches  Mehrkörpersystem  besteht  aus  einer 
Anzahl  von  Starrkörpern,  die  durch  Gelenke  (kinema¬ 
tische  Bindungen)  geometrisch,  oder  durch  Feder-/ 
Dämpferelemente  (physikalische  Bindungen)  über 
Kräfte  miteinander  gekoppelt  sind.  Starrkörper  besit¬ 
zen  eine  Geometrie  und  eine  definierte  Dichte,  über 
die  der  Schwerpunkt  des  Körpers  definiert  ist.  Nach 
dem  Satz  von  Euler  ist  die  Bewegung  der  Starrkörper 
durch  die  Bahnbewegung  (Position)  und  die  Drehbe¬ 
wegung  (Orientierung)  definiert.  Position  und  die 
Orientierung  eines  Körpers  im  Raum  können  z.B. 
durch  einen  globalen  Positionsvektor  P  und  einen 
lokalen  Rotationsvektor  R  beschrieben  werden. 

p=(p>  Pi  U 

R  =  ((p  $  i//)T 
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Ein  Körper  besitzt  daher  grundsätzlich  6  unabhängige 
Bewegungsmöglichkeiten  (Freiheitsgrade).  Die  Ge¬ 
lenke  stellen  Zwangsbedingungen  dar,  welche  die 
Freiheitsgrade  bei  den  Bewegungsmöglichkeiten  der 
Körper  einschränken  und  damit  Zwangskräfte  auf  die 
Körper  ausüben.  Sie  bestimmen  die  Kinematik  des 
Systems,  während  die  physikalischen  Bindungen  für 
das  dynamische  Verhalten  maßgeblich  sind.  Die  Si¬ 
mulation  elastischer  Fager  setzt  z.B.  voraus,  dass  die 
Körper  nicht  durch  Sperrung  von  Freiheitsgraden 
verbunden  sind,  sondern  dass  jede  Relativbewegung 
der  Körper  eine  definierte  Kraft  bewirkt.  Der  Zu¬ 
sammenhang  zwischen  Relativbewegung  und  Kraft¬ 
wirkung  kann  z.B.  durch  funktionale  Ausdrücke  oder 
Kennlinien  hergestellt  werden. 

Die  Bewegungsgleichungen  können  z.B.  beim  Ansatz 
nach  Newton-Euler  durch  das  Kräfte-  und  Momen¬ 
tengleichgewicht  formuliert  werden: 

F  =  m-ä 

M  =  [/]  •  co  +  co  x  [/]  •  cd 

Trägheitsmatrix  [7]  und  Momentvektor  M  sind  hier 
auf  den  Massenmittelpunkt  der  Körper  bezogen. 
Beim  Aufstellen  der  Bewegungsgleichungen  wird 
neben  der  „klassische  Mechanik“  auch  das  Prinzip 
von  D’Alembert  (Zwangskräfte  verrichten  keine  Ar¬ 
beit)  und  die  Erweiterung  auf  nicht  konservative  Sys¬ 
teme  (Fagrangesche  Bewegungsgleichungen)  ange¬ 
wendet.  Dadurch  werden  meist  weniger  Parameter 
und  Gleichungen  erforderlich.  Die  Fösung  des  so  er¬ 
mittelten  Gleichungssystems  kann  symbolisch  zur 
Reduktion  der  Systemordnung  als  auch  numerisch  zur 
quantitativen  Analyse  erfolgen. 

MKS  Programmsysteme  stellen  aus  dem  geometrisch 
definierten  Modell  und  den  Zwangsbedingungen 
selbständig  das  Gleichungssystem  auf  und  berechnen 
das  kinematische  und  dynamische  Systemverhalten 
im  Zeitbereich.  Als  Ergebnis  stehen  alle  Bewegungs¬ 
größen  und  innere  und  äußere  Kräfte  zur  Verfügung. 
Der  Bewegungsablauf  kann  zusätzlich  animiert  wer¬ 
den.  Zumeist  werden  im  Automobilbereich  kom¬ 
merzielle  Programmsysteme  wie  z.B.  ADAMS, 
DADS  oder  Simpack  eingesetzt. 

In  der  Fahrzeugtechnik  fand  die  Mehrkörpersimulati¬ 
on  erst  relativ  spät  Anwendung.  Hemmnisse  waren 
das  teilweise  stark  nichtlineare  Verhalten  von  Federn 
und  Dämpfern,  die  komplexe  Elastokinematik  und 
die  Schwierigkeit  bei  der  ausreichend  exakten  Model¬ 
lierung  des  Reifens  für  fahrdynamische  Untersuchun¬ 
gen.  Heute  werden  zusätzlich  zu  den  „klassischen“ 
MKS -Elementen  auch  stark  nichtlineare  Gummilager, 
dynamische  Hydrolager,  fluiddynamische  Dämpfer¬ 
modelle,  flexible  Körper  sowie  weitere  systemdyna¬ 
misch  komplexe  Komponenten  einbezogen.  Bei  der 
Modellierung  komplexer  Komponenten  kommen  so¬ 
wohl  physikalisch  begründete  als  auch  empirische 


Modelle  zum  Einsatz.  Im  ersteren  Fall  erfolgt  die 
Nachbildung  der  Komponente  in  einem  komplexen 
Ersatzmodell  aus  der  Bauteilgeometrie  und  Material¬ 
eigenschaften  (Elastizitäten,  Dämpfungsparametern 
und  Massenbelägen),  während  im  zweiten  Fall  eine 
Approximation  der  versuchstechnisch  ermittelten 
Kennlinien  durch  mathematische  Funktionen  erfolgt 
(Kennlinien-Modelle). 

Angesichts  der  hohen  Modellierungsgüte  der  Fahr¬ 
werks-  und  Gesamtfahrzeugmodelle  stellt  sich  mehr 
denn  je  die  Frage  nach  einer  geeigneten  Reifenbe¬ 
schreibung,  da  der  Reifen  wie  kein  anderes  Bauteil 
das  dynamische  Systemverhalten  mitbestimmt.  Bis¬ 
her  gibt  es  aber  kein  Reifenmodell,  das  für  alle 
Fragestellungen  gleichermaßen  gut  geeignet  ist.  Auf 
der  einen  Seite  existieren  physikalisch  orientierte 
Schwingungsmodelle  und  auf  der  anderen  Seite  empi¬ 
risch  ermittelte  Handlingmodelle.  Am  populärsten  ist 
sicher  die  „Magic  Formula“  von  Pacejka  [2].  Diese 
ist  eine  rein  mathematische  Beschreibung  des  quasi¬ 
stationären  Systemverhaltens  und  damit  für  hochdy¬ 
namische  Situationen  nur  bedingt  anwendbar.  Für 
Komfortuntersuchungen  werden  physikalisch  be¬ 
gründete  Modelle  benötigt,  die  auch  in  Bereichen 
valide  sind,  die  nicht  durch  Messungen  belegt  sind 
und  die  auch  bei  kurzwelligen  Straßenunebenheiten 
und  Einzelhindernissen  angewendet  werden  können. 
2D-Komfort-Schwingungsmodelle  (wie  z.B.  CTyre) 
können  naturgemäß  keine  Seitenkräfte  und  damit  kei¬ 
ne  Handling-  und  Sturzeinflüsse  darstellen,  während 
echte  3D-Strukturmodelle  (DTyre,  FTyre)  sehr  re¬ 
chenzeitintensiv  sind  [1]. 

15.2  Struktur-Optimierung 

Im  Bestreben  der  Fahrzeughersteller,  die  Entwick¬ 
lungszeiten  zu  verkürzen,  spielen  CA-Methoden  eine 
immer  größere  Rolle.  Ein  Schlagwort  hierbei  ist 
Upfront-CAE  -  die  zeitliche  Vorverlagerung  des  Ein¬ 
satzes  von  CA-Methoden  in  der  funktionalen  Fahr¬ 
zeugentwicklung.  Ziel  ist  es  so  früh  wie  möglich  ein 
stimmiges  Konzept  hinsichtlich  der  an  das  Fahrzeug 
gestellten  Anforderungen,  wie  zum  Beispiel  Schwin¬ 
gungskomfort  und  Akustik,  zu  erarbeiten.  Über¬ 
proportional  höhere  Kosten  bei  späteren  Konzept¬ 
änderungen  sowie  Untersuchungen  an  teueren  Hard¬ 
wareprototypen  werden  dabei  vermieden.  In  der 
Fahrzeugentwicklung  wird  in  diesem  Rahmen  schon 
seit  langem  die  Strukturberechnung  und  -Optimierung 
mittels  Finite-Elemente-Methoden  (FEM)  eingesetzt. 
Dies  beginnt  in  den  allerersten  Konzeptentwick¬ 
lungsphasen  mit  der  Erarbeitung  und  Bewertung  von 
Konzeptvarianten,  bis  hin  zur  detaillierten  Auslegung 
und  Optimierung  in  der  Serienentwicklung  eines  neu¬ 
en  Fahrzeugs. 

Ein  Finite-Elemente-Modell  für  die  Konzeptphase  ist 
dabei  anders  zu  gestalten  als  ein  Modell  für  die  Se¬ 
rienentwicklung.  So  wird  eine  detaillierte  CAD- 
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Geometrie  nicht  benötigt  und  steht  auch  oft  nicht  zur 
Verfügung.  Als  Grundsatz  für  die  Modellerstellung 
gilt  „so  grob  wie  möglich,  so  fein  wie  nötig“.  Unter 
diesen  Gesichtspunkten  fehlender  Geometriedefiniti¬ 
onen,  größtmöglicher  Flexibilität  in  der  Modellerstel¬ 
lung  und  der  Forderung  nach  Aussagen  in  kurzer 
Zeit,  ist  die  Simulation  mit  Balken-Schalen-FE- 
Modellen  eine  sehr  effektive  Methode. 

Das  Ziel  dieser  Konzeptuntersuchungen  und  Optimie¬ 
rungen  besteht  darin,  globale  Aussagen  für  das  Ge¬ 
samtfahrzeugkonzept  hinsichtlich  vorgegebener  funk¬ 
tionaler  Anforderungen,  insbesondere  Schwingungs¬ 
komfort  und  Akustik,  zu  treffen.  Dabei  sind  mögliche 
Gewichts-  und  Bauraumpotenziale  auszuweisen  so¬ 
wie  vergleichende  Konzeptvariantenuntersuchungen 
durchzuführen. 

Bei  der  Verwendung  von  Balken-Schalen-FE-Mo- 
dellen  erfolgt  die  Aufteilung  des  Fahrzeugs  nicht  to¬ 
pologisch,  sondern  nach  funktionalen  Gesichtspunk¬ 
ten  in  Träger  und  Blechfelder.  Ein  Beispiel  für  eine 
typische  Trägerstruktur  ist  der  Schweller  oder  auch 
die  B-Säule,  während  das  Fahrzeugdach  eine  typische 
Blechfeldstruktur  darstellt.  Die  Blechfeldstrukturen 
werden  mit  Schalenelementen  (2D-Elementen)  und 
die  Trägerstrukturen  mit  Balkenelementen  (1D- 
Elementen)  modelliert.  In  Bild  15-3  ist  ein  Beispiel 
eines  derartigen  Balken-Schalen-FE-Modells  einer 
Fahrzeugkarosserie  dargestellt. 


Bild  15-3:  Balken-Schalen-Konzeptmodell  eines 
Cabriolets 

Diese  Balken-Schalen-FE-Modelle  sind  besonders 
gut  geeignet  für  so  genannte  Parameteroptimierun¬ 
gen.  Bei  den  Schalenelementen  kann  nur  die  Scha¬ 
lendicke  geändert  werden,  während  bei  den  Balken¬ 
elementen  die  gesamte  Balkenquerschnittsgeometrie 
über  verschiedene  beschreibende  Parameter,  wie 
Querschnittshöhe,  -breite  und  zugeordnete  Wandstär¬ 
ken  verändert  werden  kann.  Dies  ermöglicht  einen 
wesentlich  größeren  Freiraum  für  Design-Änderun¬ 
gen  bei  den  Optimierungszyklen.  Beispiele  für  para- 
metrisierte  Balkenquerschnitte  sind  in  Bild  15-4  dar¬ 
gestellt.  Bislang  war  eine  Abstrahierung  auf  äquiva¬ 


lente  Standardquerschnitte,  wie  zum  Beispiel  Recht¬ 
eck-,  U-  oder  F-Profile,  für  die  Parametrisierung 
notwendig  Der  äquivalente  Standard-Rechteckquer¬ 
schnitt  abstrahiert  den  realen  Balkenquerschnitt  auf 
ein  Rechteckprofil  mit  den  Beschreibungsparametern 
Höhe  H,  Breite  B  und  Blechstärke  r.  Seit  kurzem  ist 
auch  eine  geometrisch  nahezu  exakte  parametrisierte 
Darstellung  des  Balkenquerschnitts  im  Rahmen  der 
Konzeptoptimierung  möglich. 


Äquivalenter  Geometrie  des  Realer 
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Bild  15-4:  Parametrisierte  Balkenquerschnitts¬ 
geometrie 

Aus  dem  Wirkkettenverständnis  sämtlicher  kundenre¬ 
levanten  funktionalen  Anforderungen  sind  in  der 
Fahrzeugentwicklung  verschiedene  Standardlastfälle 
für  die  Auslegung  und  Optimierung  der  Karosserie¬ 
struktur  definiert  worden.  Wichtig  hinsichtlich 
Schwingungskomfort  und  Akustik  sind  dabei  z.B.  die 
globalen  statischen  und  dynamischen  Karosseriestei¬ 
figkeiten  (Kapitel  4.4).  Hohe  statische  Steifigkeit  geht 
einher  mit  hoher  dynamischer  Steifigkeit,  die  sich  in 
der  Qualität  des  Schwingungsverhaltens  eines  Fahr¬ 
zeugs  ausdrückt.  Die  globale  statische  Steifigkeit 
wird  anhand  globaler  Biege-  und  Torsionslastfälle 
bewertet,  wobei  für  die  globale  dynamische  Steifig¬ 
keit  die  Frequenzlage  und  Schwingform  der  ersten 
globalen  Eigenschwingformen  überwacht  werden. 
Ein  weiterer  Fastfall  ist  der  Schwingungspegel  am 
Fenkrad.  Um  durch  Resonanzen  mit  der  Struktur  ver¬ 
ursachte  Fenkradvibrationen  gering  zu  halten,  darf 
die  Amplitude  der  Fenkradschwingung  einen  Maxi¬ 
malwert  nicht  überschreiten. 

Zu  diesen  Steifigkeits-  und  NVH-Fastfällen  können 
außerdem  noch  verschiedene  linearisierte  Crash- 
Fastfälle  wie  Frontalcrash  und  Überschlag  hinzu¬ 
kommen.  Ziel  hierbei  ist  die  Struktur  so  auszulegen, 
dass  im  Crash-Fall  die  Fahrgastzelle  erhalten  bleibt 
und  nicht  kollabiert. 

Als  nächster  Schritt  erfolgt  der  Aufbau  des  Optimie¬ 
rungsmodells  (Design-Modell).  Dazu  sind  zunächst 
die  Design- Variablen,  Querschnittshöhe,  -Breite  und 
Blechwandstärken  zu  erstellen  und  mit  den  Balken¬ 
querschnitten  zu  verknüpfen.  Die  minimalen  und  ma¬ 
ximalen  Grenzen  der  Design- Variablen  können  dabei 
von  produktionstechnischen  Randbedingungen,  wie 
zum  Beispiel  minimale  fertigbare  Blechwandstärke, 
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Bauraumrestriktionen  oder  funktionalen  Anforderun¬ 
gen  bestimmt  sein.  Für  die  komplette  Karosserie¬ 
struktur  ergeben  sich  somit  ca.  1500  Design- 
Variablen.  Hinzu  kommen  einige  tausend  Konstruk¬ 
tionsrandbedingungen  (Design-Constraints)  wie  z.B. 
Wandstärkensprünge  oder  Höhen-Seitenverhältnisse. 
Die  Design- Variablen  zusammen  mit  den  Design- 
Constraints  stellen  den  ,  Freiraum’  für  die  Optimie¬ 
rung  da. 

Das  Optimierungsmodell  wird  dann  vervollständigt 
mit  der  Formulierung  einer  Zielfunktion.  Ziel  der  Op¬ 
timierung  ist  im  Design-Freiraum  der  Struktur  eine 
optimale  Struktur  zu  ermitteln,  die  bei  minimalem 
Gewicht  alle  Anforderungen  und  Nebenbedingungen 
der  einzelnen  Lastfälle  erfüllt.  Dafür  wird  die  Viel¬ 
zahl  von  Nebenbedingungen  zu  Kategorien  zusam¬ 
mengefasst  und  zusammen  mit  den  einzelnen  Teilzie¬ 
len  mit  Straftermen  und  Gewichtungsfaktoren  in  einer 
globalen  Zielfunktion  dargestellt. 

Eine  Parameteroptimierung  eines  solchen  Karosserie¬ 
konzeptmodells  liefert  in  der  Regel  eine  große  Da¬ 
tenmenge  an  Ergebnissen  die  über  spezifische  grafi¬ 
sche  Auswerteverfahren  zu  interpretieren  sind.  In 
Bild  15-5  sind  zum  Beispiel  die  Veränderungen  der 
Balkenquerschnittsdimensionen  nach  der  Optimie¬ 
rung  relativ  zum  Ausgangszustand  dargestellt.  Deut¬ 
lich  ersichtlich  sind  die  Bereiche  der  Struktur,  in  de¬ 
nen  in  Form  von  Bauraum-  oder  Blech  Wandstärken¬ 
veränderungen  entweder  Verstärkungsmaßnahmen 
notwendig  sind  oder  Materialeinsparpotenziale  ge¬ 
funden  wurden. 


Bild  15-5:  Veränderung  der  Balkenquerschnitts¬ 
dimensionen  nach  Optimierung,  relativ  zum  Aus¬ 
gangszustand 


Mit  diesen  Ergebnissen  der  Auslegungsoptimierung  - 
in  der  Regel  zusammen  mit  Potenzial-  und  Schwach¬ 
stellenanalysen  an  verschiedener  Konzeptvarianten  - 
ist  es  nun  möglich,  in  einer  sehr  frühe  Entwicklungs¬ 
phase  wichtige  Empfehlungen  für  das  Gesamtfahr¬ 
zeugkonzept  hinsichtlich  einer  Vielfalt  von  kundenre¬ 
levanten  Anforderungen  zu  geben. 


15.3  Akustik-Berechnung  (FEM) 

Die  verfügbare  Performance  bei  der  Computersimula¬ 
tion  erlaubt  mittlerweile  Netzfeinheiten  und  -detail- 
lierungen  in  den  Finite-Elemente  Modellen,  die  -  im 
Prinzip  -  eine  Simulation  der  akustischen  Eigenschaf¬ 
ten  von  Fahrzeugen  bis  in  den  Frequenzbereich  von 
mehreren  100  Hz  möglich  macht.  Die  CAD-gestützte 
Erstellung  entsprechender  FE-Modelle  von  Karosse¬ 
riestruktur,  Interieur,  Achsen  und  Antrieb  mit  insge¬ 
samt  weit  über  zwei  Millionen  Elementen  ist  heute 
Stand  der  Technik  (Bild  15-6). 


Bild  15-6:  FE-Feinmodell  eines  PKW  der  Mittel¬ 
klasse  mit  2  Mio.  Elementen,  Struktur  und  Innen¬ 
raumkavität 


Dies  ist  die  positive  Sicht.  Richtig  ist  nämlich  auch, 
dass  es  nach  heutigem  Stand  nicht  möglich  ist,  das 
vom  Passagier  wahrgenommene  Fahrgeräusch  auch 
nur  für  einen  Betriebsfall  im  Detail  quantitativ  zu 
prognostizieren.  Mehr  noch,  dies  ist  nicht  einmal  ein 
sinnvolles  Ziel.  Eine  solche  Prognose  würde  nämlich 
einen  Reifegrad  in  der  Konstruktion,  sowie  in  der 
Beschreibung  von  frequenz-  und  temperaturabhän¬ 
gigen  Materialeigenschaften  erfordern,  die  eine 
sinnvolle  Verwendung  der  Ergebnisse  für  den  Ent¬ 
wicklungsprozess  schon  rein  zeitlich  unmöglich 
macht. 

Die  Kunst  eines  möglichst  wertschöpfenden  Entwick¬ 
lungsprozesses  besteht  also  in  der  geeigneten  Balance 
zwischen  Detaillierung  und  Prognosegüte  auf  der 
einen  Seite,  und  Robustheit  gegenüber  Unschärfen, 
Zeit-  und  Kostenzwängen  im  Projekt  auf  der  anderen 
Seite. 

Ausgehend  von  den  mathematischen  Grundlagen  der 
Akustiksimulation  soll  dieser  Sachverhalt  anhand 
einiger  in  diese  Simulation  eingehenden  physikali¬ 
schen  Größen  vertieft  werden.  Im  Anschluss  werden 
dann  Möglichkeiten  eines  Auslegungsprozesses  auf¬ 
gezeigt,  welche  zwar  nicht  mehr  direkt  das  Fahr¬ 
geräusch  adressieren,  aber  im  Sinne  einer  robusten 
und  vor  allem  rechtzeitigen  Aussage  geeignet  sind, 
die  Fahrzeugauslegung  zielführend  zu  unterstüt¬ 
zen. 
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15.3.1  Modale  Gleichungen 

Die  Simulation  des  Körperschallanteiles  im  Fahrzeu¬ 
ginnengeräusch  basiert  auf  gekoppelten  Finite- 
Elemenete-Modellen  der  Karosseriestruktur  und  des 
Luftfluids.  Die  Bewegungsgleichungen  dafür  lassen 
sich  wie  folgt  darstellen: 
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Dabei  bezeichnet  M  die  Masse-Matrix,  B  die 
Dämpfungs-  und  K  die  Steifigkeits-Matrix  von  Struk¬ 
tur  (Index  s)  und  Fluid  (Index  f).  U  sind  die  Ver¬ 
schiebungen  der  Struktur,  P  die  Drücke  im  Fluid  und 
F  die  äußeren  Kräfte  an  den  Krafteinleitungsstellen, 
also  den  Anbindungspunkten  von  Antrieb  und  Fahr¬ 
werk.  Die  Matrix  C  beinhaltet  im  wesentlichen  die 
Flächennormalen  an  den  Kontaktflächen  zwischen 
Struktur  und  Fluid  und  koppelt  damit  deren  Bewe¬ 
gungsgleichungen. 


ks,F-Q>s,F  Ks,f  *&s,f 
c  =  Of  C  Os 


mit  den  modalen  Massen-,  Steifigkeits-  und  Dämp¬ 
fungsmatrizen  von  Struktur  und  Fluid.  Erstere  beide 
sind  reine  Diagonalmatrizen.  Für  den  i-ten  Platz  auf 
der  Diagonale  gilt: 


2 


mit  der  i-ten  Eigenfrequenz  von  Struktur,  bzw.  Fluid. 

Einsetzen  ergibt  schließlich  die  gekoppelten  modalen 
Gleichungen: 
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Die  Gleichung  lässt  sich  symmetrisieren.  Dazu  diffe¬ 
renziert  man  die  erste  Zeile  nach  der  Zeit,  multipli¬ 
ziert  die  zweite  Zeile  mit  -1  und  substituiert  die  Zeit¬ 
ableitung  der  Verschiebung  U  durch  die  Geschwin¬ 
digkeit  V.  Nach  anschließender  Fourier- 
Transformation  (symbolisiert  durch  A)  erhält  man  die 
Bestimmungsgleichungen  für  die  akustisch  relevanten 
Größen  P  im  Fluid  und  Körperschallschnelle  V  in  der 
Struktur,  co  ist  die  Kreisfrequenz. 
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Damit  sind  die  Bewegungsgleichungen  des  gekoppel¬ 
ten  Fluid- Struktur  Systemes  in  ein  algebraisches 
Gleichungssystem  mit  symmetrischen  und  positiv 
definiten  Systemmatrizen  überführt.  Die  Größe  dieser 
Matrizen,  und  damit  die  Anzahl  der  zu  lösenden 
Gleichungen  ist  gleich  der  Anzahl  der  im  FE-Modell 
abgebildeten  Freiheitsgrade.  Dies  sind  oft  einige  Mil¬ 
lionen.  Darüber  hinaus  ist  dieses  Gleichungssystem 
für  jede  Frequenz  gesondert  zu  lösen. 

Durch  Transformation  in  den  Modalraum  lässt  sich 
die  Dimension  dieses  Gleichungssystemes  auf  die 
Anzahl  der  im  untersuchten  Frequenzbereich  relevan¬ 
ten  Eigenschwingungsformen  reduzieren.  Dazu  sind 
die  Eigenvektoren  und  Of  der  ungekoppelten 

Subsysteme  Struktur  und  Fluid  zu  ermitteln.  Die 
Transformation  in  den  Modalraum  erfolgt  dann  über 
die  Gleichungen 


Ü  ~  <&£  ü  >  P  ~  <&s  P  ’  f  ~  Os  F 

T 

mSf  =  *&s,F  Ms,f  *&s,f 

T 

bs,F  =  <S>s,F  Bs,f  &s,f 


Basis  dieser  Gleichungen  sind  die  Eigenvektoren 
und  .  Als  Beispiel  dazu  sind  in  Bild  15-7  jeweils 

die  ersten  vier  globalen  Eigenschwingungsformen 
von  Struktur  und  Fluid  dargestellt. 


1  Biege  eigenform 


t.  To-fsionsfiigenform 


2.  Siegeeigenform 


Vordewagemranken 


Bild  15-7:  Globale  Eigenschwingungsformen  von 
Struktur  (oben)  und  Fluid  (unten) 


15.3  Akustik-Berechnung  (FEM) 
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Der  Rechenaufwand  hängt  in  erster  Linie  von  der 
Anzahl  der  für  die  Lösung  der  modalen  Gleichung 
benötigten  Eigenmoden  ab.  Als  grobe  Faustregel 
kann  gelten,  dass  zur  Simulation  des  akustischen 
Verhaltens  bei  einer  Frequenz  /0  mindestens  alle  Ei¬ 
genmoden  mit  einer  Eigenfrequenz  kleiner  als  1 ,5/0 
ermittelt  und  zur  Lösung  der  gekoppelten  modalen 
Gleichung  herangezogen  werden.  Wie  in  Bild  15-8  zu 
sehen,  überwiegen  für  Frequenzen  unterhalb  von 
1000  Hz  die  Eigenmoden  der  Struktur  bei  weitem. 
Als  Beispiel  ist  hier  das  gekoppelte  Schwingungs¬ 
verhalten  bei  200  Hz  unter  Berücksichtigung  von 
400  Struktur-  und  20  Fluidmoden  angegeben. 


5U00- 


Bild  15-8:  Eigenmoden  in  Abhängigkeit  von  der  Fre¬ 
quenz  mit  einem  Beispiel  für  eine  gekoppelte  Simula¬ 
tion  bei  200  Hz 

Vorstehend  wurde  hergeleitet,  wie  bei  gegebenen  Ei¬ 
genschaften  des  Systems  -  aus  Sicht  der  Simulation 
also  bei  gegebenen  Massen-,  Dämpfungs-  und  Steifig¬ 
keitsmatrizen  -  das  Akustik-  und  Schwingungsver¬ 
halten  berechnet  werden  kann.  Für  homogene,  metalli¬ 
sche  Werkstoffe  ist  die  Breitstellung  ausreichend  ge¬ 
nauer  Daten  hierzu  auch  kein  größeres  Problem.  Die 
Herausforderungen  stecken  in  Inhomogenitäten  an  den 
Verbindungsstellen,  sowie  in  den  temperatur-  und  fre¬ 
quenzabhängigen  Materialeigenschaften  nichtmetalli¬ 


scher  Werkstoffe.  Hiervon  ist  neben  der  Steifigkeit  vor 
allem  die  Dämpfung  betroffen  (Kapitel  6.2.6). 


•  • 

15.3.2  Äquivalente  abgestrahlte 
Schallleistung 

Die  Entstehung  des  Fahrzeuginnengeräusches  (ge¬ 
nauer,  dessen  Körperschallanteil)  lässt  sich  als  eine 
Kette  von  Wirkungsmechanismen  darstellen.  Jeder 
dieser  Mechanismen  kann  in  Blockschaltbildform  - 
mit  jeweils  definierter  Zuordnung  von  Eingangs-  und 
Ausgangskenngrößen  der  verschiedenen  Blöcke  - 
dargestellt  werden.  Diese  Übertragungspfade  sind 
von  primärer  Bedeutung  hinsichtlich  der  Ge¬ 
räuschübertragung  im  Frequenzbereich  der  dominan¬ 
ten  Motorordnungen,  also  etwa  bis  500  Hz. 


Bild  15-9:  Wirkkette  für  den  Körperschallanteil  am 
F  ahrzeuginnengeräusch 

Mit  Bezug  auf  diese  Übertragungskette  kann  für  die 
im  Fahrzeuginnenraum  erzeugten  Schalldrücke  p(a>) 
geschrieben  werden: 

p(co)  =  H_(co)  ■  c(co)  •  G(co)  •  F{(jo)  , 
mit 

F(&>):  Vektor  der  frequenzabhängigen  Betriebserre¬ 
gerkräfte  im  Antriebsstrangsystem,  wie  zum  Beispiel 
Trägheits-  und  Gaswechselkräfte,  sowie  der  aus  Stra¬ 
ßen-  und  Unwuchterregung  resultierenden  Kräfte, 

G(co) :  Frequenzgangmatrix  des  Übertragung s Verhal¬ 
tens  zwischen  den  Betriebserregerkräften  des  An¬ 
triebsstranges  und  den  Auslenkungen  der  antriebs¬ 
strangseitigen  Koppelpunkte  zur  Karosserie, 

c(co) :  frequenzabhängige  Diagonalmatrix  der  dyna¬ 
mischen  Steifigkeiten  der  Lager  zwischen  Antriebs¬ 
strang  und  Karosserie, 

:  Frequenzgangmatrix  des  mechanisch-akusti¬ 
schen  Übertragung s Verhaltens  zwischen  den  Karosse¬ 
rie-Anregungskräften  und  den  Schalldrücken  an  defi¬ 
nierten  Positionen  im  Fahrzeuginnenraum. 
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Alle  oben  aufgeführten  Größen  sind  komplex,  das 
heißt  phasenbehaftet.  Es  kann  gezeigt  werden,  dass 
schon  durch  Änderung  der  wechselseitigen  Phasen¬ 
beziehungen  zwischen  den  einzelnen  Primärkräften  in 
Motor  und  Fahrwerk  Streubänder  im  Schalldruck  von 
etwa  6  dB  entstehen. 

Zur  Reduzierung  dieser  Komplexität  ist  es  gängige 
Praxis,  je  nach  Phase  und  Reifegrad  des  Projektes  ge¬ 
eignete  Hilfs-  und  Komponentenziele  zur  Berücksich¬ 
tigung  der  akustischen  Anforderungen  zur  definieren. 
Diese  orientieren  sich  an  der  oben  beschriebenen 
Wirkkette.  Als  Beispiel  soll  hier  auf  das  Hilfsziel 
Äquivalente  abgestrahlte  Schallleistung4  näher  einge¬ 
gangen  werden.  Dabei  geht  man  von  der  maschinen¬ 
akustischen  Gleichung  aus  (Kapitel  6.1.4): 

P{co)  =  Z0-S-vi-a(ct>) 

vn  ist  die  Normalkomponente  der  Schwingschnelle 
auf  der  Oberfläche  der  Karosseriestruktur,  S  ist  die 
Ab  strahlfläche  und  o  ist  der  Abstrahlgrad.  Werden 
die  Beiträge  der  einzelnen  infinitesimalen  abstrahlen¬ 
den  Blechflächen  addiert,  d.h.  die  Auslöschungsef¬ 
fekte  aufgrund  der  Phasenbeziehungen  vernachläs¬ 
sigt,  so  kann  der  Abstrahlgrad  er  =1  gesetzt  werden. 
Damit  definiert  sich  die  äquivalente  abgestrahlte 
Schallleistung  (ERP  für  equivalent  radiated  power) 
einer  Karosseriesubstruktur  als: 

ERP(co)  =  Z0  •  co2  •  ||ww(eu)|2  ds 

Der  Vorteil  dieses  Vorgehens  besteht  zum  einen  da¬ 
rin,  dass  zur  Simulation  der  ERP-Kurven  akustische 
Effekte  wie  Absorption  und  Auslöschung  nicht  be¬ 
trachtet  werden,  und  damit  auch  keine  Annahmen 
darüber  getroffen  werden  müssen.  Weiterhin  ermög¬ 
licht  dies  eine  robuste  Mittelung  der  gewichteten 
Schwinggeschwindigkeiten  über  sinnvolle  Karosse¬ 
riebaugruppen  (Bild  15-10). 


Als  Ergebnis  dieses  Verfahrens  erhält  man  einen 
,ERP-Fingerprint‘  der  Karosseriestruktur  (Bild 
15-11):  Im  unteren  Diagramm  sind  hier  im  Frequenz¬ 
bereich  bis  400  Hz  die  äquivalenten  abgestrahlten 
Schallleistungen  für  einzelne  Karosseriebaugruppen, 
sowie  deren  Summe  dargestellt.  Dies  liefert  einen 
schnellen  Überblick  über  eventuell  kritische  Blech¬ 
feldschwingen.  Das  obere  Diagramm  zeigt  die  zuge¬ 
hörigen  berechneten  Schalldruckkurven  für  die  Pas¬ 
sagierplätze,  sowie  die  Eigenfrequenzen  der  Fluid¬ 
moden  des  Fahrzeuginnenraumes. 


Bild  15-11:  Beispiel  ERP-, Fingerprint4  einer  Karos¬ 
seriestruktur 


15.4  Boundary  Element  Methode 
(BEM) 

Bei  der  Boudary  Element  Methode  (BEM)  wird  im 
Gegensatz  zur  „klassischen44  Finite  Elemente  Metho¬ 
de  im  Rahmen  der  Diskretisierung  des  Problems  nicht 
der  ganze  Betrachtungsraum  sondern  nur  die  Oberflä¬ 
che  modelliert  (Randwertverfahren).  Dadurch  kann 
das  Problem  um  eine  Ordnung  (3D  — ►  2D)  reduziert 
und  die  Anzahl  der  Freiheitsgrade  drastisch  einge¬ 
schränkt  werden.  Damit  wird  eine  einfache  Aufberei¬ 
tung  des  Rechennetzes  aus  den  CAD-Daten  möglich 
[4]. 


Heckfensterrahmen 


Stirnwand 


Heckklappe 


Frontscheibe 


Dach 


Bild  15-10:  Beispiel  für  die  Unterteilung  der  Karos¬ 
serie  in  Baugruppen  zur  ERP- Simulation 


FEM  (Gebietsverfahren)  BEM  (Randwertverfahren) 


EirueNast 


Bild  15-12:  Gegenüberstellung  Diskretisierung 
FEM/BEM 

In  einem  ersten  Schritt  werden  nun  die  Randwerte  auf 
der  Oberfläche  (z.B.  Schalldruck  oder  Schallschnelle) 
bestimmt.  In  einem  zweiten  Schritt,  ausgehend  von 
diesen  Randbedingungen,  die  Schallausbreitung  im 
umgebenden  Fluid  auf  Basis  der  bekannten  Funda- 
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287 


mentallösung  für  die  Schallausbreitung  ermittelt  (Ka¬ 
pitel  5.2.2).  Mit  dieser  Methode  ist  die  Modellierung 
und  Berechnung  unendlich  ausgedehnter  Gebiete 
leicht  möglich,  was  speziell  bei  der  Schallausbreitung 
im  Luftraum  außerhalb  der  Fahrgastzelle  relevant  ist. 
Dementsprechend  wird  BEM  vorteilhaft  insbesondere 
für  die  Berechnung  des  Vorbeifahrgeräusches  einge¬ 
setzt. 

Das  Vorbeifahrgeräusch  eines  Fahrzeugs  setzt  sich 
aus  den  Beiträgen  Antriebsgeräusch  und  Reifen- 
Fahrbahngeräusch  zusammen  (Kapitel  14).  Abwei¬ 
chend  von  der  rein  energetischen  Summation  der 
Einzelbeiträge  erzeugen  die  weitgehend  kohärenten 
Teilquellen  von  Ansaug-  und  Schalldämpfermündung 
ein  interferierendes  Schallfeld  und  damit  Welligkei¬ 
ten  im  Vorbeifahrpegel.  Mit  der  Transferpfadanalyse 
(Kapitel  16.4)  kann  der  Schalldruckpegel  pM  am  Mik¬ 
rofonpunkt  phasenrichtig  aus  dem  Schalldruck  der 
Ansaug-  und  der  Schalldämpfermündung  ( pAM „  pSM) 
auf  addiert  werden: 

Pm  ~  PaM  '  TFAM  +  P SM  ■  TFsm 

Die  Transferfunktionen  von  der  Ansaug-  und  Schall¬ 
dämpfermündung  zum  Mikrofonpunkt  (TFAM,  TFSM ) 
können  sowohl  experimentell  als  auch  rechnerisch 
ermittelt  werden.  An  Stelle  der  Messung  der  FRFs 
bietet  sich  deren  Berechnung  mittels  der  Boundry 
Element  Methode  an.  Das  BEM-Netz  wird  dabei 
abhängig  von  dem  interessierenden  Frequenzbereich 
in  angemessener  Feinheit  über  die  Fahrzeugstruk¬ 
tur  gelegt.  Die  Oberfläche  wird  als  schallhart  defi¬ 
niert. 

Beispiele  von  so  berechneten  FRFs  in  Bild  15-14  für 
die  Abstrahlung  von  der  Schalldämpfermündung  zu 
ausgewählten  Aufpunkten  an  den  Vorbeifahrmikro¬ 
fonen  zeigen,  dass  der  charakteristische  Frequenz¬ 
gang  in  Rechnung  und  Messung  gut  harmoniert.  Die 
verbleibenden  Abweichungen  in  der  Frequenzlage 
von  Maxima  und  Minima  sind  unerheblich  für  die 
Bewertung  durch  Terz-  oder  Summenpegel.  Bei  der 
Synthese  des  Schallpegelverlaufs  für  einen  Mikro- 


Bild  15-13:  BEM-Netz  der  Oberflächenstruktur  eines 
Fahrzeugs 


Bild  15-14:  Gemessene  und  mit  BEM  berechnete 
FRFs  im  Vergleich 


fonpunkt  im  Außengeräusch  ergeben  sich  damit  aus¬ 
reichend  gute  Übereinstimmungen  zwischen  Messung 
und  Rechnung  (Bild  15-16).  Die  verbleibenden  Ab¬ 
weichungen  sind  nicht  zuletzt  durch  die  Restreflekti- 
vität  des  Semi-Freifeldraums  verursacht. 


In  der  Fahrzeugkonzeption  kann  diese  Berechnungs¬ 
methode  zur  Bewertung  von  Vorbeifahrpegeln,  etwa 
zur  Beurteilung  des  Einflusses  der  (virtuellen)  Fahr- 


Bild  15-15:  Mit  BEM  berechnete  horizontale  Schall¬ 
feldverteilung  des  Mündungsgeräusches  bei  30  Hz 
(oben)  und  115  Hz  (unten) 
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zeugform,  herangezogen  werden,  wie  die  beobachte¬ 
ten  Beugungseffekte  für  eine  simulierte  Schalldämp¬ 
fermündungsquelle  nahe  legen  (Bild  15-15).  Offen¬ 
sichtlich  werden  bei  dieser  Schallfeldverteilung  die 
Kavitäten  im  Bereich  des  Motorraumes  und  des  Un¬ 
terbodens  an  die  Umgebung  des  Fahrzeugs  angebun¬ 
den.  Die  Modellierungsgüte  ist  an  die  jeweilige  Fra¬ 
gestellung  anzupassen. 


Bild  15-16:  Vergleich  Vorbeifahrgeräusch  gemessen 
und  berechnet  mit  BEM 


Bei  der  rechnerischen  Auslegung  von  Abgasanlagen 
stehen  die  durch  den  Gaswechsel  angeregten  Mün¬ 
dungsschalle  zur  Verfügung.  Mit  der  hier  vorgestell¬ 
ten  Methode  kann  deren  Abstrahlung  ebenfalls  rech¬ 
nerisch  weiterverfolgt  werden.  Im  Berechnungsmo¬ 
dell  besteht  zusätzlich  die  weitergehende  Möglichkeit 
der  Anbindung  von  Oberflächenschnellen  der  an¬ 
grenzenden  Strukturen  (gewonnen  etwa  aus  FEM- 
Analysen)  als  zusätzliche  Luftschallquellen.  Dem 
entsprechend  kann  der  berechnete  Schalldruck  an  den 
begrenzenden  Flächen  des  Fahrzeuginnenraumes 
auch  zur  Bewertung  von  Einflüssen  des  Außenschalls 
auf  die  Innenraumakustik  dienen. 

15.5  Statistische  Energieanalyse 
(SEA) 

Die  statistische  Energieanalyse  (SEA)  ist,  wie  der 
Name  schon  ausdrückt,  ein  statistisches  Berech¬ 
nung  s  verfahren,  bei  dem  mit  Erwartungswerten  der 
gespeicherten  Energien  und  ausgetauschten  Leistun¬ 
gen  zwischen  gekoppelten  schwingungsfähigen  Sys¬ 
temen  gerechnet  wird.  Auf  eine  exakte  Berechnung 
der  Schallwellen  nach  Betrag  und  Phase,  wie  es  in  der 
wellentheoretischen  Akustik  üblich  ist,  wird  dabei 
verzichtet.  Es  wird  vielmehr  in  einem  betrachteten 
Subsystem  von  einer  gewissen  Anzahl  von  Schwin¬ 
gung  smoden  mit  der  Modendichte  AN/Af  und  dem 


gesamten  Energieinhalt  W  ausgegangen.  Um  eine 
gewisse  statistische  Sicherheit  zu  erreichen,  muss 
über  möglichst  viele  Moden  gemittelt  werden.  Des¬ 
halb  wird  über  einen  ausreichen  großen  Frequenzbe¬ 
reich,  gewöhnlich  Terzen  oder  Oktaven,  gemittelt.  In 
diesem  Bereich  sollte  eine  Mindestanzahl  von  Reso¬ 
nanzfrequenzen  liegen.  Diese  Mindestanzahl  wird  in 
der  Literatur  mit  sehr  unterschiedlichen  Werten  ange¬ 
geben,  ein  guter  Anhaltspunkt  ist  dabei  der  Wert  von 
AN  =  6  Moden  pro  Frequenzband  Af.  Bei  dreidimen¬ 
sionalen  Kavitäten,  wie  sie  bei  der  Luftschallübertra¬ 
gung  vorliegen,  steigt  die  Modendichte  nach  Kapitel 
5.6  quadratisch  mit  der  Frequenz  an: 


AN 

V 


«4  K'V ' 


c 


3 

0 


Berücksichtigt  man  den  Zusammenhang  zwischen 
Bandmittenfrequenz  fn  und  Bandbreite  Af  bei  Terz- 
und  Oktavfiltern: 

Af  =  0,71-  fm  Oktavfilter 

Af  =  0,11-  fm  Terzfilter 

so  kann  die  notwendige  Mittenfrequenz  bestimmt 
werden,  damit  innerhalb  der  Filterbandbreite  die  er¬ 
forderliche  Anzahl  von  Resonanzen  enthalten  sind 
(Bild  15-17).  Die  Anwendbarkeit  der  SEA  im  Luft¬ 
schallbereich  ist  damit  auf  Frequenzbereiche  be¬ 
schränkt,  die  oberhalb  dieser  Filterfrequenzen  liegen. 
Je  kleiner  die  Kavitäten  werden,  und  je  schmalbandi- 
ger  die  Betrachtung  erfolgen  soll,  umso  höher  liegt 
diese  Frequenzgrenze.  In  der  Praxis  kann  von  einem 
nutzbaren  Frequenzbereich  oberhalb  von  400-500  Hz 
ausgegangen  werden. 


Bild  15-17:  Bandmittenfrequenz  von  Terz-  und 
Oktavfrequenzbändem  die  gerade  6  Moden  enthalten 

Ein  einfaches  Beispiel  für  die  Methodik  der  SEA  ist 
in  Bild  15-18  dargestellt.  Die  Elemente  1  und  3  stel¬ 
len  dort  Sende-  und  Empfangsraum  dar,  während 
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Element  2  das  Schwingsystem  Trennwand  repräsen¬ 
tiert.  Der  Verlustfaktor  rjn  beschreibt  in  diesem  Fall 
die  resonante  Anregbarkeit  von  Biegewellen  in  der 
Trennwand,  rj23  die  Abstrahlung  der  Trennwand  in 
den  Empfangsraum  und  rjl3  die  nicht  resonante  Über¬ 
tragung  von  Raum  1  nach  Raum  3  infolge  des  Mas¬ 
sengesetzes. 


7ll  =7l+7l2+7l3 
bzw.  allgemein 

>!:  ■  -  ~  Z  '!■ 

i*i 


Bild  15-18:  SEA-Modell  einer  Trennwand  (2)  zwi¬ 
schen  Senderaum  (1)  und  Empfangsraum  (3) 


Wenn  z.B.  in  dem  Körper  1  keine  Leistung  Px  mehr 
zugeführt  wird,  erfolgt  infolge  der  dissipativen  Ver¬ 
luste  eine  exponentielle  Abnahme  der  gespeicherten 
Energie.  Die  Verluste  werden  durch  den  Dämpfungs¬ 
verlustfaktor  rjx  charakterisiert: 

dW 

dt 

Bei  zwei  gekoppelten  multimodalen  Schwingsyste¬ 
men  mit  der  Modenanzahl  N1  und  N2  wird  Schwin¬ 
gungsenergie  von  Mode  zu  Mode  übertragen.  Die  pro 
Schwingung  ausgetauschte  Energie  ist  proportional 
zur  Differenz  der  modalen  Energien: 


R 


l<-»2 


oj-N1-N2 


rwL_wp 

vV  äT 


Die  zwischen  System  1  und  2  übertragene  Leistung 
kann  damit  auch  in  Abhängigkeit  von  Kopplungsver¬ 
lustfaktoren  ;//2und  //^formuliert  werden: 

P^2  =  ^ '  (7i2  ‘  W,  —  tj21  ■  kV2 )  =  Pl2  —  P2l 

Die  Energiebilanz  für  das  System  von  Bild  15-18 
ergibt  sich  damit  wie  folgt: 


'p' 

7  ii  -72i  -7si 

R 

=  CO- 

— 7l2  722  “732 

• 

w2 

c, 

_— 7l3  “723  733 

Wo 

(P)  =  o3\r1\(W) 


Mit  dem  SEA-Modell  kann  nun  der  Energiepegel 
jedes  Subsystems  berechnet  werden,  wenn  man  die 
Verlustfaktoren  und  die  eingespeiste  Leistung  kennt. 
Als  Ergebnis  der  Berechnungen  erhält  man  die  ge¬ 
samte  Schwingungsenergie  im  betrachteten  Fre¬ 
quenzband.  Daraus  können  dann  die  interessierenden 
dynamischen  Größen  berechnet  werden.  Mathema¬ 
tisch  erhält  man  den  Energiepegel,  indem  man  For¬ 
mel  15-1  nach  der  Energie  auflöst: 

W=+W'U) 

CO 

Für  den  Zusammenhang  zwischen  dem  räumlich  ge¬ 
mittelten  Effektiv  wert  des  Schalldrucks  und  der 
Energie  in  einer  Kavität  mit  dem  Volumen  V  und  der 
charakteristischen  Impedanz  z0  erhält  man: 


Unter  der  Annahme,  dass  nur  im  Senderaum  1  Ener¬ 
gie  eingespeist  wird  (P2  =  P3  =  0),  errechnet  sich  z.B. 
aus  obigen  Gleichungen  für  den  Isolationspegel  zwi¬ 
schen  System  1  und  System  3: 

W 

L,  =  101g-—/-  =  201g 

vv  o 


L  (%  +  Jhi  +  %3  )■(%+%  i )  +  (%  +  % i )  rin 

rhi{rh+rhi+rin)  +  rhi%i 

Da  die  SEA  die  Subsysteme  durch  stark  vereinfachte 
Annahmen  beschreibt  und  dadurch  die  Anzahl  der 
Parameter  vergleichsweise  klein  bleibt,  kann  man 
z.B.  Sensitivitätsanalysen  recht  einfach  durchführen 
und  so  dominante  Schallpfade  erkennen  bzw.  das 
Systemverhalten  optimieren.  Dies  ist  ein  erheblicher 
Vorteil  gegenüber  anderen  numerischen  Berech¬ 
nungsmethoden.  Allerdings  ist  die  Unsicherheit  der 
Analyse  infolge  der  statistischen  Berechnungsgrund¬ 
lagen  auch  größer. 

SEA-Anwendungen  können  in  zwei  Gruppen  unter¬ 
teilt  werden: 


mit  dem  Gesamtverlustfaktor  von  System  1: 


(15-1) 


Experimentelle  Untersuchung  von  Modellen 

Die  Anwendung  der  analytischen  SEA  für  Prog¬ 
nose-Rechnungen 
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Die  experimentelle  Untersuchung  erfolgt  hauptsäch¬ 
lich  mit  der  sogenannten  Power  Injection  Method 
(PIM).  Die  PIM  erlaubt  die  Bestimmung  von  SEA- 
Parametern  in  situ,  ohne  dass  man  die  Subsysteme 
entkoppeln  muss,  sowie  die  Bestimmung  von  Trans¬ 
missionswegen.  Grundlage  dieser  experimentellen 
SEA  ist  wie  bei  der  analytischen  SEA  die  Leistungs¬ 
bilanz  der  einzelnen  Subsysteme.  Zur  Bestimmung 
der  Eingangsparameter  wird  Formel  15-1  jedoch  nach 
der  Matrix  [r\]  aufgelöst,  die  die  Dämpfungs-  und 
Kopplungsverlustfaktoren  beinhaltet: 

CO 


Daraus  ergibt  sich  bereits,  dass  bei  der  Messung  die 
Energien  jedes  Subsystems  sowie  die  in  jedes  Sub¬ 
system  eingeleitete  Leistung  bestimmt  werden  müs¬ 
sen.  Normalerweise  kann  die  Energie  nicht  direkt 
gemessen  werden,  sondern  es  wird  die  mittlere  Be¬ 
schleunigung  einer  Struktur  oder  der  mittlere  Schall¬ 
druckpegel  einer  Kavität  erfasst. 

Für  die  üblichen  Dämmungs-,  Dämpfungs-  und  Ab¬ 
sorptionsmaterialien  existieren  bereits  diverse  Ma¬ 
terialmodelle,  mit  deren  Hilfe  die  Dämpfungsver¬ 
lustfaktoren  auf  physikalische  Grunddaten,  wie  z.B. 
Dichte,  Elastizitätsmodul,  Strömungswiderstand, 
usw.  zurückgeführt  werden  können  (Tabelle  15-1). 
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Die  Kopplungen  zwischen  Lufträumen  (Kavitäten) 
und  Festkörpern  sind  detailliert  untersucht  und  die 
Kopplungsverlustfaktoren  gut  bekannt.  Kopplungen 
zwischen  steifen  Strukturen  sind  dagegen  wesentlich 
problematischer,  da  der  Energieaustausch  bei  höheren 
Frequenzen  stark  von  geometrischen  Details  abhängt. 
Daher  ist  die  Anwendung  der  SEA  im  Körperschall¬ 
bereich  auch  heute  noch  eher  Gegenstand  von  metho¬ 
dischen  Grundsatzuntersuchungen,  während  sie  sich 
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Tabelle  15-1:  Materialdaten  und  Materialmodelle  für 
SEA  Anwendungen  im  Luftschallbereich 


Bild  15-19:  SEA-Ergebnisse  und  Modell  zur  Luft- 
schall-Übertragung  vom  Motorraum  in  den  Innen¬ 
raum 

bei  der  Auslegung  von  Luftschallpfaden  im  Automo¬ 
bil  mittlerweile  als  Standardmethode  etabliert  hat. 
Mit  AutoSEA  steht  dafür  auch  ein  kommerzielles 
SW- Werkzeug  zur  Verfügung. 

Die  rechnerische  Betrachtung  der  Luftschallpfade  mit 
Hilfe  der  SEA  ermöglicht  eine  optimale  Auslegung 
der  Schallisolierung  bereits  in  frühen  Projektphasen. 
Mit  entsprechendem  Modellierungsaufwand  kann 
dabei  z.B.  der  Energiefluss  von  den  Eingangskavitä¬ 
ten  im  Motorraum  bis  zum  Innenraum  der  Fahrgast¬ 
zelle  genau  beschrieben  werden.  Im  SEA-Modell  las¬ 
sen  sich  so  neben  dem  Beitrag  einzelner  Stirn  wand¬ 
bereiche  zum  Dämmung sverhalten  auch  Leckageein¬ 
flüsse  analysieren  (Bild  15-19).  Mittlerweile  sind 
auch  Optimierungsalgorithmen  verfügbar,  mit  deren 
Hilfe  eine  weitgehend  automatisierte  Verteilung  von 
Dämmungs-  und  Dämpfungsbelägen  hinsichtlich  Di¬ 
cke  und  Anordnung  erfolgen  kann.  Insbesondere  in 
Verbindung  mit  spritzbaren  Belägen  kann  so  eine 
gewichtsoptimale  Lösung  realisiert  werden. 

15.6  Aeroakustik-Berechnung  (CAA) 

Die  Computational  Aeroacoustic  (CAA)  benutzt  Me¬ 
thoden,  die  das  Umströmungsgeräusch  mit  unter- 


15.6  Aeroakustik-Berechnung  (CAA) 


291 


schiedlichen,  vereinfachenden  Annahmen  numerisch 
berechnen.  Diese  Methoden  liefern  Aussagen  über 
Wechseldrücke  auf  den  Oberflächen  und  Fluktuati¬ 
onen  in  der  Umströmung,  die  signifikante  akustische 
Quellen  darstellen.  Dabei  bilden  CFD- Verfahren 
(Computational  Fluid  Dynamics)  die  Grundlage  für 
die  numerische  Aeroakustik.  Diese  berechnen  Strö¬ 
mungsgeschwindigkeiten  v(/)  und  Drücke  p{t).  In  der 
Akustik  sind  nur  die  Schwankungsanteile  v  und  p 
von  Interesse  (Schallschnelle  und  Schalldruck).  Eine 
besondere  Herausforderung  stellt  dabei  die  Tatsache 
dar,  dass  sich  die  typischen  aerodynamischen  und  die 
akustischen  Skalen  stark  unterscheiden.  Man  beden¬ 
ke,  dass  die  Schmerzgrenze  des  menschlichen  Ohres 
bei  Wechseldrücken  von  100  Pa  liegt.  Bei  einem 
Fahrzeug  tritt  im  Vergleich  schon  ein  statischer  Un¬ 
terdrück  von  mindestens  1 15  Pa  bei  50  km/h  auf.  Au¬ 
ßerdem  ist  das  Hörempfinden  logarithmisch,  so  dass 
schon  kleine  Wechseldrücke  von  0.02  Pa  als  Sprech¬ 
lautstärke  wahrgenommen  werden  (Kapitel  7).  Bei 
einer  Strömung  mit  Ma  =  0,1  wird  deshalb  z.B.  eine 
Genauigkeit  der  Strömungsberechnung  von  min¬ 
destens  10"9  benötigt.  Wird  diese  nicht  erreicht,  so 
äußern  sich  die  Ungenauigkeiten  der  numerischen 
Berechnung  als  „numerischer  Lärm“.  Darüber 
hinaus  gibt  es  große  Disparitäten  der  Längenskalen 
zwischen  dem  Wirbelmaßstab  der  Strömung  und 
der  Wellenlänge  des  durch  die  Wirbel  erzeugten 
Schalls. 

Die  Methoden  zur  Simulation  des  Umströmungsge¬ 
räusches  können  in  die  Direkte  Numerische  Simula¬ 
tion  (DNS)  und  in  hybride  Ansätze  unterteilt  werden. 
Die  Basis  für  die  direkte  numerische  Simulation  bil¬ 
den  die  kompressiblen  Navier-Stokes-Gleichungen 
(N-S).  Es  werden  sowohl  die  strömungsrelevanten  als 
auch  die  akustischen  Skalen  betrachtet  und  auch  die 
turbulenten  Skalen  direkt  berechnet.  Jedoch  erzwingt 
die  Auflösung  der  einzelnen  Turbulenzen  ein  sehr 
feines  Gitter,  so  dass  dies  eigentlich  nur  in  der 
Forschung  unter  Zuhilfenahme  von  Supercomputem 
und  bei  kleinen  Reynolds-Zahlen  sinnvoll  möglich 
ist. 

In  der  Praxis  haben  sich  daher  aus  Aufwandsgründen 
hybride  Lösungen  auf  Basis  der  Reynolds-Gleichungen 
durchgesetzt.  Diese  beschreiben  zeitlich  gemittelte  Zu¬ 
standsgrößen  des  Strömungsfeldes.  Deshalb  ist  in  die¬ 
sem  Fall  ein  separates  Turbulenzmodell  nötig,  um  das 
Gleichungssystem  zu  schließen.  Als  Mittelweg  gilt  die 
Large  Eddy  Simulation,  die  zumindest  die  großen  Wir¬ 
bel  direkt  numerisch  berechnet  und  erst  die  kleinen 
Skalen  über  ein  Turbulenzmodell  simuliert. 

Ein  neuerer  Ansatz  ist  die  numerische  Simulation  auf 
der  Grundlage  der  Lattice-Boltzmann-Methode 
(LBM).  Dabei  wird  das  Fluid  nicht  als  Kontinuum 
sondern  als  Teilchen- Strom  aufgefasst  und  dessen 
Mikrodynamik  berechnet  (Kapitel  15.6.2). 


Aus  den  direkt  berechneten  oder  modellierten  Druck- 
und  Geschwindigkeits-Fluktuationen,  die  als  Quell¬ 
terme  dienen,  wird  dann  mit  Hilfe  von  linearisierten 
Euler-Gleichungen  (LEE)  oder  mit  Hilfe  akustischer 
Analogien  das  Schallfeld  im  Fernfeld  bestimmt.  Zu 
den  gängigen  akustischen  Analogien  zählen  u.a.  die 
Verfahren  nach  Lighthill  (Kapitel  5.2.1),  Ffowcs  Wil¬ 
liams  &  Hawkings. 


Newtonsches 

Fluid 


, - 


Kontinuumsmechanik 


kompressible 
N-S-Gleichungen 
•  " — 

Ma<0,3 


V 

Teilchenmechanik 

V._ 

A 


Lattice-  Boltzmann- 


Bild  15-20:  Schema  der  aerodynamischen  und  aero- 
akustischen  Berechnung 


Für  die  CAA  stehen  mittlerweile  auch  einige  leis¬ 
tungsfähige  kommerzielle  Softwaretools  zur  Verfü¬ 
gung.  Insbesondere  das  auf  der  LBM  (Kapitel  15.6.2) 
basierende  PowerFLOW®  von  ©Exa  Corporation,  ist 
für  die  Akustiksimulation  geeignet.  Fluent®  von 
©ANSYS,  Ine.  und  STAR-CD™  von  CD-adapco 
Group  nutzen  dagegen  N-S -basierte  Algorithmen  der 
Kontinuumsmechanik.  Weiterhin  werden  an  Univer¬ 
sitäten  und  Forschungsstätten  verschiedene  Pro¬ 
gramme  entwickelt,  die  zum  einen  das  Ziel  haben, 
den  durch  die  Strömung  direkt  verursachten  Lärm, 
zum  anderen  aber  auch  den  durch  strömungsangereg¬ 
te  Oberflächen  erzeugten  Lärm  zu  erfassen. 


15.6.1  Navier-Stokes-Gleichungen 

Bei  den  bisher  gängigsten  Methoden  der  numerischen 
Strömungsmechanik,  wie  z.B.  die  Direkte  Numeri¬ 
sche  Simulation,  die  Reynolds  Averaged  Navier- 
Stokes-Gleichungen  oder  die  Large  Eddy  Simulation, 
wird  das  newtonsche  Fluid  als  Kontinuum  betrachtet. 
Aus  den  Erhaltungsgleichungen  für  Masse,  Impuls 
und  Energie  entsteht  ein  System  gekoppelter,  nichtli¬ 
nearer,  partieller  Differentialgleichungen  2.  Ordnung, 
das  auch  die  Turbulenz  und  die  hydrodynamische 
Grenzschicht  enthält. 
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Die  Navier-Stokes-Gleichungen  sind  die  Grundglei¬ 
chungen  der  Strömungsmechanik.  Sie  setzen  sich  aus 
der  Kontinuitätsgleichung  und  dem  Impulserhal¬ 
tungssatz  zusammen  und  beschreiben  die  Strömungs¬ 
geschwindigkeit  und  die  Druckverteilung  in  newton- 
schen  Flüssigkeiten.  Die  allgemeinste  Form  der  Na¬ 
vier-Stokes-Gleichungen  lautet: 


du 


+ 


p(w-V)w  =  —Vp  +  rjÄu  +  ('A,  +  r})V{y-u')  +  f 


Der  Vektor  /beschreibt  die  Volumenkraftdichte  (z.B. 
die  Gravitation  oder  die  Corioliskraft),  p  den  physika¬ 
lischen  Druck.  Die  Viskositäten  2  und  rj  werden  als 
bekannt  vorausgesetzt.  Ist  die  Dichte  entlang  der 
Teilchenbahnen  konstant,  ergeben  sich  mit  der  Kon¬ 
tinuitätsgleichung 

^-  +  V\pu)  =  0 

dt  '  ’ 


•  Finite-Elemente-Methode  (FEM). 

•  Lattice-Boltzmann-Methode  (LBM) 

•  Randelementmethode  (boundary  element 
method,  BEM) 

Da  es  sich  bei  allen  oben  angeführten  Methoden  um 
numerische  Näherungs  verfahren  handelt,  ist  eine  Va¬ 
lidierung  mit  dem  Versuch  in  jedem  Fall  erforderlich. 
Bei  turbulenten  Strömungen  treten  für  die  numerische 
Strömungssimulation  noch  einige  Fragen  auf.  Entwe¬ 
der  werden  sehr  feine  Rechengitter  wie  bei  der  DNS, 
oder  zusätzliche  Turbulenzmodelle  notwendig,  wel¬ 
che  neben  den  numerischen  Fehlern  auch  noch  Mo¬ 
dellierung  sfehler  verursachen.  Hier  besitzt  die  LBM 
prinzipielle  Vorteile  gegenüber  FE-basierten  Metho¬ 
den. 

15.6.2  Lattice-Bolzmann-Methode 


die  Divergenzfreiheit  des  Geschwindigkeitsfeldes,  die 
ein  inkompressibles  Fluid  charakterisiert  und  damit 
die  so  genannten  inkompressiblen  Navier-Stokes- 
Gleichungen.  Ist  die  Dichte  des  Fluids  konstant,  so 
vereinfachen  sich  die  Gleichungen  zu 


du 

dt 


+ 


( u-V)u  =  -Vp  +  uAu  +  f 


Dabei  wird  der  diffuse  Impulstransport  mit  der  kine¬ 
matischen  Viskosität  v  =  rj/p  beschrieben.  Diese  Glei¬ 
chungen  sind  für  Luftströmungen  weit  unterhalb  der 
Schallgeschwindigkeit  (Ma  <  0,3)  gültig. 


Eine  weitere  Vereinfachung  sind  die  Euler-Glei¬ 
chungen.  Sie  enthalten  keine  Turbulenz  und  keine 
Grenzschicht,  da  die  Reibung  nicht  mit  in  Betracht 
gezogen  wird.  Allerdings  kann  mit  diesem  Ansatz 
kein  Strömungsabriss  simuliert  werden.  Sie  eignen 
sich  daher  nur  für  Berechnungen,  in  denen  die  Grenz¬ 
schicht  keine  wesentliche  Rolle  spielt. 


Bei  der  statistischen  Turbulenz-Modellierung  wird  im 
einfachsten  Fall  die  turbulente  Strömung  als  zeitlicher 
Mittelwert  und  Varianz  der  Geschwindigkeit  v  und 
dem  Druck  p  gebildet.  Wird  diese  Annahme  in 
die  Navier-Stokes-Gleichungen  eingefügt  und  zeit¬ 
lich  gemittelt,  erhält  man  die  Reynolds  Averaged  Na¬ 
vier-Stokes-Gleichungen  (RANS)  mit  dem  Rey¬ 
nolds’ sehen  Spannungstensor.  Es  kann  anhand  der 
verwendeten  Turbulenzgrößen  zwischen  Null-,  Ein- 
und  Zweigleichungsmodellen  unterschieden  werden. 
Das  k-e-Turbulenzmodell  z.B.  ist  ein  Zweiglei¬ 
chungsmodell,  das  mit  zwei  partiellen  Differential¬ 
gleichungen  die  Entwicklung  der  turbulenten  kineti¬ 
schen  Energie  k  und  der  isotropen  Dissipationsrate  s 
beschreibt. 


Gängige  numerische  Lösungsmethoden  der  Strö¬ 
mungsmechanik  sind  die 


Die  numerische  Umsetzung  der  Boltzmann- 
Gleichung  erfolgte  Ende  der  80er  Jahre  in  Form  des 
Lattice-Boltzmann- Verfahrens.  Im  Gegensatz  zu  den 
klassischen  Verfahren  basiert  das  Lattice-Boltzmann- 
Verfahren  nicht  auf  einer  Diskretisierung  der  eine 
Strömung  beschreibenden  Navier-Stokes-Gleichung, 
sondern  auf  der  aus  der  statistischen  Mechanik  be¬ 
kannten  Boltzmann-Gleichung,  d.h.  auf  der  Berech¬ 
nung  einer  stark  vereinfachten  Teilchen-Mikrodyna¬ 
mik.  Das  Fluid  wird  als  eine  Ansammlung  kleiner 
Teilchen  beschrieben,  die  sich  analog  der  kinetischen 
Gastheorie  im  Raum  bewegen.  Es  gibt  zwei  verschie¬ 
dene  Bewegungsarten.  Eine  setzt  sich  aus  der  Advek- 
tion  und  Kollision  im  Raum  zusammen,  die  andere 
beschreibt  die  Kollision  von  Partikeln  mit  Oberflä¬ 
chen.  Die  gemittelten  Größen  weisen  somit  ein 
makroskopisches  Verhalten  auf,  das  umso  realisti¬ 
scher  wird,  je  feiner  das  Gitter  gewählt  wird.  Zusätz¬ 
lich  wird  die  Phasenraumdichte  eingeführt,  die  jedem 
Geschwindigkeitsvektor  eine  Teilchendichte  im  An¬ 
schauungsraum  zuordnet.  Die  Evolution  dieser  Pha¬ 
senraumdichte  ist  durch  die  Boltzmann  Gleichung 
gegeben.  Diese  Form  der  Beschreibung  wird  als  me- 
soskopisch  bezeichnet.  Es  werden  weder  die  einzel¬ 
nen  Teilchen  betrachtet  (mikroskopisch),  noch  wird 
das  System  auf  seine  wesentlichen  makroskopi¬ 
schen  Kenngrößen  reduziert.  In  der  Kinetik  ist  die 
Chapman-Enskog  Entwicklung  etabliert,  um  aus  der 
mesoskopischen  Boltzmann-Gleichung  die  makro¬ 
skopische  Euler-  und  Navier- Stokes  Gleichung  her¬ 
zuleiten.  Aufgrund  der  internen  Struktur  (geringer 
Speicher-  und  Rechenbedarf  je  Zelle)  eignet  sich  das 
Verfahren  u.a.  zur  Berechnung  von  Strömungen  in 
komplexen  Geometrien. 

Bild  15-21  zeigt  eine  aeroakustische  Simulation  im 
Vergleich  zur  Messung  am  Unterboden  einer  Ober¬ 
klasse-Limousine.  Die  Messungen  wurden  auf  einer 
mittleren  Höhe  zwischen  Fahrbahn  und  Unterboden 
in  einem  vorher  festgelegten  Raster  durchgeführt.  Die 
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Simulationsebene  wurde  zur  Validierung  in  die  glei¬ 
che  Ebene  gelegt.  Das  Ergebnis  von  Simulation  und 
Messung  zeigt  eine  sehr  gute  Übereinstimmung  der 
Strukturen  und  lokalen  Verteilungen  des  Schallwech- 
seldruckes  bei  40  Hz. 


Bild  15-21:  Simulation  Schallwechseldruck  im  Un¬ 
terbodenbereich  bei  40Hz  (oben)  und  Messung  im 
Windkanal  (unten) 

Weitere  wichtige  Anwendungsgebiete  für  aeroakusti- 
sche  Berechnungen  sind  die  Simulation  der  Umströ¬ 
mung  eines  Cabrioverdecks  sowie  im  Hecknachlauf 
und  die  Überströmung  von  Kavitäten.  Mit  diesen  In¬ 
formationen  können  bereits  in  der  Konzeptphase  Aus¬ 
sagen  über  den  Strömungsverlauf,  die  Wirbelbildung 
und  die  erwarteten  Schalldruckpegel  gemacht  werden. 


der  Struktur  bereits  festgelegt  ist,  z.B.  durch  Über¬ 
nahmeteile,  während  ein  weiterer  Teil  der  Struktur 
noch  zu  optimieren  ist.  In  der  Regel  sind  diese  Sys¬ 
temteile  elastisch  gekoppelt,  so  dass  das  Strukturver¬ 
halten  beider  Teile  nicht  unabhängig  voneinander 
analysiert  werden  kann.  Bei  der  FE-Modellierung 
dient  die  Substrukturtechnik  der  Eliminierung  des 
unveränderlichen  Teils  der  Strukturgleichungen  und 
verringert  so  den  numerischen  Aufwand.  Dafür  ge¬ 
eignete  Substrukturtechniken  sind  bereits  seit  Jahr¬ 
zehnten  bekannt  [6].  Oft  können  bei  der  numerischen 
Simulation  jedoch  nicht  alle  Teilsysteme  eines  kom¬ 
plexen  Systems,  wie  z.B.  die  mechanisch-akustische 
Übertragungsfunktionen  einer  Fahrzeugkarosserie, 
analytisch  ausreichend  exakt  abgebildet  werden  (Ka¬ 
pitel  15.3).  Es  besteht  deswegen  in  der  Schwingungs¬ 
technik  und  Akustik  vielfach  der  Wunsch,  numeri¬ 
sche  Modelle  von  bestimmten  Teilsystemen  mit  ge¬ 
messenen  Eigenschaften  anderer  Teilsysteme  zu  einer 
Gesamtsystemaussage  zu  kombinieren.  In  der  Litera¬ 
tur  wird  eine  derartige  Vorgehens  weise  auch  als  hyb¬ 
ride  Modellierung  bezeichnet.  Die  Kopplung  der  ge¬ 
messenen  und  simulierten  Komponenten  und  Teilsys¬ 
teme  zum  Gesamtsystem  setzt  ebenfalls  geeignete 
Substrukturtechniken  voraus. 

Substrukturtechniken  führen  das  Systemverhalten 
grundsätzlich  auf  die  fundamentalen  Bewegungsglei¬ 
chungen  der  Systeme  zurück.  Das  bedeutet,  dass  bei 
der  Kopplung  mehrerer  (Sub)Systeme  immer  folgen¬ 
de  drei  Sachverhalte  berücksichtigt  werden  müssen: 

•  Die  Übertragungseigenschaften  zwischen  den 
Ein-  und  Ausgangsgrößen  der  einzelnen  Subsys¬ 
teme; 


Es  ist  abzusehen,  dass  sich  für  die  numerische  Aero- 
akustik  ein  ähnlich  weites  Anwendungsgebiet  öffnet, 
wie  es  vor  ca.  zehn  Jahren  bei  der  Aerodynamik  ent¬ 
standen  ist.  Die  Rechnerleistung  und  der  Entwick¬ 
lungsstand  der  numerischen  Methoden  machen  es 
möglich,  die  aeroakustische  Lärmentstehung  und 
-ausbreitung  zusammen  mit  der  Strömung  zu  berech¬ 
nen.  Aus  dem  Verständnis  der  Strömungs-  und  Tur¬ 
bulenzeffekte  lassen  sich  gezielt  geräuschmindernde 
Maßnahmen  über  die  Strömungsbeeinflussung  ablei¬ 
ten.  Aeroakustische  Lärmquellen  können  somit  iden¬ 
tifiziert  und  reduziert  werden.  Auch  die  Wechselwir¬ 
kung  zwischen  Strömung  und  elastischen  Körpern 
sowie  der  dadurch  induzierte  Lärm  spielt  im  Fahr¬ 
zeugbau  eine  wichtige  Rolle. 

15.7  Substrukturtechnik1 

15.7.1  Methodische  Grundlagen 

Bei  der  Entwicklung  komplexer  mechanischer  Struk¬ 
turen  ergibt  sich  häufig  die  Situation,  dass  ein  Teil 

1  Im  angelsächsischen  Sprachraum:  dynamic  substructuring  (DS) 


•  Das  Kräftegleichgewicht  zwischen  den  Koppel¬ 
punkten  der  Subsysteme  (Equilibrium  Condi¬ 
tion); 

•  Die  Gleichheit  der  Bewegungen  an  den  Kopp¬ 
lungspunkten  zwischen  den  Subsystemen 
(Compatibility  Condition); 

Bei  der  sogenannten  direkten  Methode  (Direct 
Coupling)  werden  die  Subsysteme  im  Zeitbereich 
(Physical  Domain)  über  die  Systemparameter  Masse, 
Dämpfung  und  Steifigkeit  miteinander  gekoppelt.  Sie 
ist  damit  unmittelbar  äquivalent  zur  Standard  FE  Mo¬ 
dellierung  und  z.B.  auch  die  Grundlage  der  Modellie¬ 
rung  in  Kapitel  4.  Die  sogenannte  „Component  Mode 
Synthesis“  oder  CMS-Methode,  koppelt  dagegen  die 
verschiedenen  Subsysteme  auf  Basis  ihrer  modalen 
Parameter  (Modal  Domain).  Der  Grundgedanke  dabei 
ist,  dass  sich  das  dynamische  Verhalten  des  Gesamt¬ 
systems  durch  eine  begrenzte  Anzahl  von  Eigenmo¬ 
den  mit  geringer  Dämpfung  beschreiben  lässt.  Da  die 
Anzahl  der  Eigenmoden  z.B.  in  Luftschallkavitäten 
mit  der  Frequenz  schnell  ansteigt  (Kapitel  5.6.1),  ist 
diese  Beschreibungsmethode  in  diesem  Fall  auf  rela¬ 
tiv  niederfrequente  Phänomene  beschränkt.  Sie  findet 
insbesondere  im  Rahmen  der  experimentellen  Modal- 
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Ein  gekoppeltes  System  gemäß  Bild  15-22  lässt  sich 
allgemein  durch  folgendes  Gleichungssystem  be¬ 
schreiben: 
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_o_ 

(15-2) 


Die  Verschiebungen  im  Vektor  u  aller  Freiheitsgrade 
der  Knoten  (Koppelpunkte)  korrespondieren  mit  den 
Steifigkeiten  der  Koppelpunkte  in  der  dynamischen 
Steifigkeitsmatrix  Z.  Dies  bedeutet,  dass  Z  eine 
Blockdiagonalmatrix  darstellt,  mit  den  Submatrizen 
Z* 11  bis  Z(n)  für  die  verschiedenen  Subsysteme: 


Bild  15-22:  Abstraktes  Schema  der  mechanischen 
Kopplung  zwischen  zwei  Subsystemen 

analyse  bei  der  Beschreibung  des  dynamischen  Ver¬ 
haltens  von  Fahrzeugstrukturen  Anwendung  (Kapitel 
16.5.1).  Bei  der  Methode  des  „Frequency  Based 
Substructuring“  (FBS)  basiert  die  Kopplung  zwischen 
den  Subsystemen  hingegen  auf  errechneten  oder  ge¬ 
messenen  frequenzabhängigen  Übertragungs¬ 
funktionen  (Frequency  Response  Function  FRF). 
Diese  Methode  ist  damit  grundsätzlich  auch  für  die 
Beschreibung  des  Systemverhaltens  im  Bereich  hoher 
Modendichte  und  -dämpfung  geeignet. 
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Bild  15-23:  Äquivalente  Beschreibungsformen  des 
dynamischen  Systemverhaltens  in  unterschiedlichen 
Domänen 
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Die  Verbindungen  zwischen  den  Freiheitsgraden  der 
Subsysteme  sind  in  der  boolschen  Matrix  B  zeilen¬ 
weise  jeweils  durch  1  und  -1  repräsentiert,  ansonsten 
ist  diese  Matrix  nur  mit  Nullen  gefüllt.  Für  das  Sys¬ 
tem  in  Bild  15-22  mit  seinen  9  Freiheitsgraden  und 
drei  Verbindungen  ergibt  sich  somit  beispielhaft: 
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Über  die  Gleichung 


B -u  =  0 


wird  die  Kompatibilitätsbedingung  sichergestellt.  Das 
Kräftegleichgewicht  wird  in  Formel  15-2  nur  implizit 
durch  die  Fagrange  Multiplikatoren  2  beschrieben 
und  so  die  Subsysteme  miteinander  verbunden.  Wird 
2  aus  dieser  Gleichung  eliminiert,  so  erhält  man  den 
vereinfachten  Zusammenhang  zwischen  den  Ver¬ 
schiebungen  und  den  Anregungskräften/: 

u  =  iY-¥BT(BYBT)  '  BY\f 


Obwohl  alle  drei  Darstellungsformen  prinzipiell 
äquivalent  sind,  ist  in  der  experimentellen  Praxis  die 
Systemanalyse  und  -synthese  mittels  CMS  oder  FBS 
am  besten  geeignet.  Handelt  es  sich  um  recht  kom¬ 
plexe  Systeme  mit  großer  Modendichte,  wie  im  Fall 
der  mechanisch-akustischen  Eigenschaften  einer 
Fahrzeugkarosserie,  so  kann  nur  die  FBS  Methode 
sinnvoll  angewendet  werden,  da  die  Systemeigen¬ 
schaften  nicht  mehr  ausreichend  über  modale  Para¬ 
meter  beschrieben  werden  können.  Nachstehend  wird 
deswegen  die  sogenannte  Fagrange  Multiplier  Fre¬ 
quency  Based  Substructuring  Methode  (FM-FBS) 
näher  dargestellt. 


Die  dynamischen  Nachgiebigkeits-Matrizen  Y  können 
experimentell  recht  einfach  unter  frei-frei  Bedingung 
mittels  Shaker  oder  Impulshammer  bestimmt  werden 
(Kapitel  16.3).  Die  Gleichung  eignet  sich  damit 
grundsätzlich  gut  für  die  experimentelle  Umsetzung 
der  DS -Methode. 

Während  DS-Methoden  schon  seit  Jahrzehnten  bei 
der  numerischen  Modellierung  erfolgreich  eingesetzt 
werden,  führt  die  Anwendung  dieser  Methoden  in 
Verbindung  mit  experimentellen  Ergebnissen  jedoch 
auf  praktische  Probleme.  Diese  sind  dafür  verant¬ 
wortlich,  dass  die  Einbindung  experimenteller  Ergeb- 
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nisse  in  modellierte  Systeme  in  vielen  Fällen  nur 
teilweise  erfolgreich  war.  Diese  Probleme  lassen  sich 
in  den  folgenden  vier  Punkten  zusammenfassen: 

•  Es  wird  eine  punktförmige  Verbindung  zwi¬ 
schen  Systemen  mit  3  translatorischen  und  3  ro¬ 
tatorischen  Freiheitsgraden  angenommen.  In  der 
Praxis  ist  eine  Punktmessung  in  alle  Raumrich¬ 
tungen  allerdings  kaum  möglich.  Es  gibt  derzeit 
auch  keine  praktikablen  Möglichkeiten,  um  Tor¬ 
sionsanregungen  auf  diesen  Punkten  anzubrin¬ 
gen. 

•  Die  idealisierten  unendlich  steifen,  masselosen 
und  punktförmigen  Verbindungen  sind  kon¬ 
struktiv  nicht  darstellbar.  In  der  Realität  handelt 
es  sich  vielmehr  immer  um  flächenhafte  Ver¬ 
bindungen  mit  begrenzten  Steifigkeiten. 

•  Die  Übertragungseigenschaften  der  einzelnen 
Subsysteme  sollen  unter  frei-frei-Bedingungen 
bestimmt  werden.  Im  Experiment  wird  das  Sys¬ 
tem  aber  mit  einer  gewissen  Steifigkeit  mit  sei¬ 
ner  Umgebung  verbunden,  so  dass  damit  Fehler 
in  der  Analyse  anfallen.  Auch  die  Anbindung 
eines  Shakers,  sowie  die  Befestigung  von  Sen¬ 
soren  verusacht  unerwünschte  Rückwirkungen 
durch  Veränderung  der  Massenverhältnisse  und 
lokalen  Steifigkeiten. 

•  Die  DS-Methode  ist  generell  empfindlich  auf 
Fehler  in  den  gemessenen  FRFs.  Damit  ist  die 
Methode  anfällig  für  Messfehler  aller  Art. 

15.7.2  EMPC-Methode 

Wie  im  vorstehenden  Abschnitt  dargestellt,  ist  die 
Berücksichtigung  der  rotatorischen  Freiheitsgrade  der 
Subsystemschnittstellen  ein  entscheidender  Faktor  bei 
der  Umsetzung  der  DS-Methode.  Da  es  jedoch  nur 
eingeschränkte  Möglichkeiten  gibt,  Momentenanre- 
gungen  am  Prüfkörper  zu  applizieren  und  dessen  Ro¬ 
tation  exakt  zu  vermessen,  fehlen  bisher  meist  erfolg¬ 
reiche  Anwendungen. 

Mit  der  nachstehend  beschriebenen  Methode  wird 
eine  Ersatzlösung  vorgestellt,  welche  auf  die  direkte 
Messung  der  rotatorische  Freiheitsgrade  verzichtet 
[5].  Die  Fösung  stützt  sich  auf  den  Ansatz  der  FEM, 
welcher  eine  komplexe  Struktur  in  einfache  kleine 
Elemente  unterteilt.  Die  Anzahlt  der  definierten  Frei¬ 
heitsgrade  an  den  Substrukturschnittstellen  beschrei¬ 
ben  die  möglichen  Schnittstellenbewegungen.  Je 
mehr  Freiheitsgrade  eine  Schnittstellte  definieren, 
umso  komplexere  Verformungen  können  zwischen 
den  Subsystemen  dargestellt  werden.  Dabei  müssen 
nicht  unbedingt  rotatorische  Freiheitsgrade  definiert 
werden.  Beispielsweise  werden  Volumenelemente 
normalerweise  nur  mit  translatorischen  Freiheitsgra¬ 
den  beschrieben.  Eine  gewisse  Anzahl  von  Knoten 
mit  ausschließlich  translatorischen  Freiheitsgraden 
reicht  ebenfalls  aus,  um  auch  rotatorische  Bewegun¬ 
gen  implizit  zu  beschreiben. 


Dieser  Ansatz  kann  nun  auch  auf  die  experimentelle 
Substrukturierung  übertragen  werden.  Durch  mehrere 
translatorische  Übertragungsfunktionen  auf  mehrere 
Schnittstellenknoten  können  die  lokalen  Verformun¬ 
gen  ebenfalls  hinreichend  angenähert  werden.  Rotato¬ 
rische  Informationen  werden  auf  dieser  Weise  indi¬ 
rekt  erworben.  Dieser  Ansatz  wird  in  der  Fiteratur  als 
Equivalent  Multi  Point  Connection  Methode  (EMPC) 
beschrieben. 

Wird  die  EMPC  Methode  verfolgt,  so  muss  der  An¬ 
wender  sein  Augenmerk  auf  die  Schnittstellendefini¬ 
tion  legen.  Wenn  z.B.  Schnittstellen  zwischen  be¬ 
nachbarten  Subsystemen  an  Stellen  mit  großen  Stei¬ 
figkeitsübergängen  definiert  werden,  so  wird  die 
Schnittstelle  über  lokale  Starrkörperbewegungen 
meist  ausreichend  beschrieben.  Da  das  weichere  Sub¬ 
system  den  Bewegungen  des  steiferen  Systems  folgt, 
wird  die  Schnittstellenbeschreibung  vereinfacht,  was 
auch  zu  Vereinfachungen  bei  der  Analyse,  reduzier¬ 
tem  Messaufwand  und  besseren  Kopplungsergebnis¬ 
sen  führt. 

Als  Anwendungsbeispiel  zeigt  Bild  15-24  die  Hinter¬ 
achse  eines  PKWs.  Heutzutage  kann  der  Antrieb¬ 
strang  mit  seinen  Komponenten  wie  Hinterachsträger 
(HAT),  Getriebe  (HAG)  und  auch  die  Hinterach-  und 
Hinterachsgetriebe-Gummilager  (HAT-F  &  HAG-F) 
ziemlich  genau  modelliert  werden  (Kapitel  15.1).  Ihre 
Auswirkung  im  Gesamtfahrzeug  lässt  sich  allerdings 
modelltechnisch  nicht  befriedigend  bestimmen,  da 
sich  die  Karosserie  mit  ihrem  Übertragungsverhalten 
bis  zum  Fahrerohr  nur  sehr  eingeschränkt  rechnerisch 
darstellen  lässt  (Kapitel  15.3).  Deshalb  bietet  sich 
hier  die  hybride  DS-Methode  an,  um  das  numerische 
Antrieb  Strangmodell  mit  gemessenen  Übertragungs¬ 
funktionen  der  Karosserie  zu  kombinieren.  Die  Kopp¬ 
lung  ermöglicht  damit  eine  Optimierung  der  Antrieb¬ 
strangkomponenten  in  der  realen  Gesamtfahrzeug¬ 
umgebung. 


Bild  15-24:  Hinterer  Antriebstrang  eines  Fahrzeugs 
mit  Standardantrieb  und  Fagerpunkten 
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Das  numerische  Model  des  Hinterachsträgers  kann 
z.B.  in  NASTRAN  analysiert  werden.  Mittels  frei-frei 
Randbedingungen  werden  die  Eigenfrequenzen  und 
Moden  im  Frequenzbereich  von  0  bis  1500  Hz  be¬ 
rechnet.  Mittels  einer  modalen  Synthese  wird  daraus 
die  Übertragungsfunktion  des  Hinterachsträgers  er¬ 
stellt.  Für  das  akustische  Verhalten  stellen  die  Gum¬ 
milager  wegen  ihrer  Isolationseigenschaften  die 
wichtigsten  Komponenten  des  Antrieb  Strangs  dar.  Da 
die  experimentelle  Bestimmung  der  Fagereigenschaf- 
ten  unter  Berücksichtigung  des  rotatorischen  Verhal¬ 
tens,  wie  vorstehend  bereits  diskutiert,  nur  beschränkt 
möglich  ist,  können  diese  durch  entsprechend  detail¬ 
lierte  Fagermodelle  ermittelt  werden  (Bild  15-25). 
Diese  berücksichtigen  sowohl  die  großen  Fagerver- 
formungen  beim  Einbau  als  auch  die  frequenzabhän¬ 
gige  Steifigkeit  und  Dämpfung  gemessener  Material¬ 
proben.  Die  äußeren  und  inneren  Aluminiumteile 
werden  dabei  als  Starrkörper  modelliert,  da  ihre  Stei¬ 
figkeit  gegenüber  dem  Gummi  um  ein  Vielfaches 
größer  ist. 


HMHi  t  bnliirn  HAT-L  HAG4,  Vorne 


Bild  15-25:  Fagermodelle 


Für  die  experimentelle  Bestimmung  der  mechanisch¬ 
akustischen  Übertragungsfunktionen  der  Karosserie 
wird  am  Fahrzeug  die  Hinterachse  ausgebaut  und 
diese  über  Fuftfedern  abgestützt.  Die  vier  Karosserie- 
Anbindungspunkte  werden  dann  jeweils  an  drei  Kno¬ 
ten  in  drei  Raumrichtungen  mittels  Shaker  vermes¬ 
sen. 

In  Bild  15-26  sind  typische  Ergebnisse  aus  dieser 
Vorgehens  weise  dargestellt.  Dabei  zeigt  die  schwarze 
Kurve  die  Validationsmessung,  während  rot  und  blau 
SPC  und  EMPC  Kopplungsresultate  dokumentieren. 
Die  Single  Point  Connection  (SPC)  Kopplung  be¬ 
rücksichtigt  nur  die  drei  translatorischen  Freiheits¬ 
grade  pro  Schnittstelle.  Das  Bild  zeigt,  dass  die 
EMPC-Methode  im  Vergleich  zum  SPC- Ansatz  zu 
deutlich  besseren  Prognosen  für  das  Gesamtsystem 
führt.  Eine  zu  geringe  Dämpfung  und  meist  zu  hohe 
simulierte  Eigenfrequenzen  ab  400  Hz  sind  allerdings 
auch  für  die  EMPC-Methode  erkennbar.  Fetzterer 
Mangel  hat  maßgeblich  mit  der  Definition  der  starren 
Schnittstellen  zu  tun.  Die  reale  Struktur  besitzt  offen¬ 
sichtlich  mehr  Flexibilität,  was  zur  Senkung  der 
Eigenfrequenzen  beiträgt.  Die  Resultate  zeigen  deut¬ 
lich,  dass  die  rotatorischen  Freiheitsgrade  nicht  ver¬ 
nachlässigt  werden  dürfen.  Deshalb  führt  im  Allge¬ 
meinen  nur  die  EMPC-Methode  zu  befriedigenden 
Kopplungsresultaten. 


200  400  600  B00  1ÜG0 

Frequency  [Hz) 

Bild  15-26:  Gesamtsystem-Übertragungsfunktion, 
gemessen  sowie  mittels  SPC-  und  EMPC-Methode 
prognostiziert 

15.7.3  Identifikation  von 
Betriebsanregungen 

Neben  dem  Problem  der  exakten  Modellierung  des 
akustischen  Übertragungsverhaltens  von  Fahrzeugen 
stellt  die  exakte  Bestimmung  der  Betriebsanregungen 
häufig  ein  zentrales  Problem  dar.  So  können  bei¬ 
spielsweise  die  teilweise  sehr  komplexen  Anregungs¬ 
mechanismen  von  Motor  und  Getriebe  (Kapitel  9)  im 
Regelfall  nur  unvollkommen  modelliert  werden.  Dies 
lässt  sich  auch  dadurch  zu  erklären,  dass  die  Anregun¬ 
gen  häufig  stark  nichtlinear  vom  jeweiligen  Betriebs¬ 
zustand  des  Systems  abhängen,  die  Modelle  jedoch 
empfindlich  auf  Änderungen  in  der  Anregung  reagie¬ 
ren.  Um  trotzdem  zu  brauchbaren  Vorhersagen  für  das 
akustische  Verhalten  eines  Fahrzeugs  zu  kommen, 
wird  in  diesem  Abschnitt  eine  analytisch-experimen¬ 
telle  Vorgehensweise  vorgestellt,  mit  deren  Hilfe  auch 
nicht  direkt  messbare  Anregungen  ausreichend  be¬ 
rücksichtigt  werden  können  [5]. 

Die  Methode  wird  anhand  der  Darstellungen  von  Bild 
15-27  am  Beispiel  der  Betriebsanregungen,  welche 
beim  Abrollen  der  Zahnräder  in  einer  Getriebestufe 


I  1 


(a)  (b)  (c) 


Bild  15-27:  Indirekte  Identifikation  von  Betriebs¬ 
anregungen  (gb  =  gearbox,  bw  =  body  in  white) 
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auftreten,  erläutert.  Diese  Betriebsanregungen  könn¬ 
ten  nur  mit  erheblichem  experimentellem  Aufwand 
direkt  gemessen  werden.  Im  vorliegenden  Abschnitt 
wird  daher  gezeigt,  wie  an  Stelle  der  anregenden 
Zahnkraft  auch  äquivalente  Kräfte  an  messtechnisch 
besser  zugänglichen  Stellen  des  Systems  genutzt 
werden  können,  um  daraus  das  Gesamtsystemverhal¬ 
ten  abzuleiten.  Dabei  werden  Linearität  und  Zeitinva¬ 
rianz  der  Systeme  unterstellt. 

Im  Teilbild  (a)  ist  das  Gesamtsystem  vereinfachend  als 
Zweimassenschwinger  dargestellt.  Der  Ansatz  kann 
jedoch  analog  auf  eine  beliebige  Anzahl  von  Freiheits¬ 
graden  erweitert  werden.  Mathematisch  wird  das  dar¬ 
gestellte  System  mit  Hilfe  der  dynamischen  Steifigkei¬ 
ten  an  den  Schnittstellen  wie  folgt  beschrieben: 
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Dabei  ist  anzumerken,  dass  sich  das  Getriebe  am 
Prüfstand  infolge  der  geänderten  Lagerung  anders 

verhält  wie  im  realen  Fahrzeug,  was  mit  u\  berück¬ 
sichtigt  wurde.  Durch  Elimination  der  Zahnbewegung 
u\  in  Formel  (15-5)  ergibt  sich: 

fl=Z2x-  u\  =  Z21  •  Zj",1  •  fg  (15-6) 

Da/  der  wirksamen  Kraft  wie  in  Formel  (15-4)  äqui¬ 
valent  ist,  kann  diese  ersetzt  werden.  Nach  weiteren 
Umformungen  [3]  ergibt  sich  schließlich  folgende 
vereinfachte  Darstellung: 
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(15-7) 


uv  u2  und  u3  sind  die  Verschiebungen  an  der  Zahn¬ 
flanke,  am  Lagerflansch  des  Getriebes,  sowie  die 
Verschiebung  der  Lageraufnahme  an  der  Karosserie. 
Die  Verzahnungskraft  des  Getriebes  /  stellt  die  Be¬ 
triebsanregung  des  Systems  dar.  Zn  bis  Z33  sind  die 
dynamischen  Steifigkeiten  von  Getriebe  und  Karosse¬ 
rie  an  den  jeweiligen  Koppelpunkten.  Alle  Größen 
sind  frequenzabhängig  und  basieren  auf  der  bekann¬ 
ten  Bewegungsgleichung  für  lineare  Systeme  im  Fre¬ 
quenzbereich: 

Z  — co  •  A4  +  j  co  •  C  +  K 


Y22  und  F32  stellen  dabei  die  dynamischen  Nachgie¬ 
bigkeiten  des  Gesamtfahrzeugs  dar  und  können  auf 
einfache  Weise  mittels  Impulshammer  oder  Shaker 
gemessen  werden.  Die  Ableitung  macht  deutlich,  dass 
die  Verzahnungskräfte  nicht  explizit  bekannt  sein 
müssen,  um  das  Betriebs  verhalten  des  Gesamtsys¬ 
tems  zu  bestimmen.  Das  gleiche  Resultat  wird  enthal¬ 
ten,  wenn  die  Betriebskräfte  des  Getriebes  auf  einem 
Komponentenprüfstand  mit  starrer  Einspannung  ge¬ 
messen  werden  und  daraus  gemäß  Formel  15-7  das 
Betriebs  verhalten  des  ursprünglichen  Systems  be¬ 
rechnet  wird. 


Soll  an  Stelle  der  Verschiebung  an  der  Karosserie 
z.B.  der  Luftschallpegel  am  Fahrerohr  betrachtet 
werden,  so  kann  u3  durch  den  Schalldruck  p  ersetzt 
werden.  Entsprechend  ändern  sich  dann  auch  die  Ein¬ 
heiten  der  Matrixelemente  Z23,  Z32  und  Z33.  Die  Ge¬ 
samtsteifigkeit  in  den  Verbindungspunkten  errechnet 
sich  aus  der  Summation  der  Steifigkeiten  der  gekop¬ 
pelten  Teilsysteme: 

7 (tot)  _  7(gb)  ,  y(bw) 

^22  ~  ^22  ~1~  ^22 

wobei  die  Indices  (gb)  auf  das  Getriebe  und  (bw)  auf 
die  Karosserie  hinweisen.  Eliminiert  man  in  Formel 
15-3  die  Bewegung  ux  in,  so  erhält  man  die  reduzierte 
Form: 
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Gemäß  Bild  (b)  wird  das  Getriebe  nun  auf  einem 
Komponentenprüfstand  mit  ideal  steifen  Lagerpunk¬ 
ten  befestigt.  Die  Bewegungsgleichung  dieses  Sys¬ 
tems  kann  wie  folgt  formuliert  werden: 


In  der  Praxis  erweist  es  sich  jedoch  zumeist  als 
schwierig,  die  Lageraufnahmen  des  Komponenten¬ 
prüfstandes  im  relevanten  Frequenzbereich  von  100- 
1000Hz  in  Übereinstimmung  mit  obiger  Annahme 
ausreichend  steif  zu  gestalten.  Deshalb  wird  nachste¬ 
hend  eine  Erweiterung  der  Methode  vorgenommen, 
um  zusätzlich  auch  noch  die  Prüfstandseinflüsse  zu 
eliminieren. 
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Bild  15-28:  Schematische  Darstellung  der  Situation 
beim  Prüfstand  mit  elastischen  Lagerpunkten 
(tb  =  testbench) 

Der  Ansatz  ist  in  Bild  15-28  schematisch  dargestellt. 
Wie  bisher  wird  das  Getriebe  von  den  unbekannten 
Verzahnungskräften  /  angeregt.  Für  diese  Konfigura- 
tion  können  folgende  Gleichungssysteme  formuliert 
werden: 
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Die  Sensorkraft  fs  ist  dabei  die  Kraft  im  Lagerpunkt 
des  Getriebes  auf  dem  Prüfstand.  Wird  ux  aus  der  ers¬ 
ten  Gleichung  eliminiert,  so  enthalt  man: 

(z<f  -z21z-'z12)M2  =  I  -z2lz;lfs 

Verknüpft  man  dieses  Resultat  mit  Formel  (15-6),  so 
ergibt  sich  nach  einigen  Umformungen: 


nrK  oear  orcer  an&nuae  an  RAD  n  ImefTare  n  y  d  recöm 


Bild  15-29:  Synthese  des  Schalldrucks  am  Fahrerohr 
(unten)  aus  den  am  Prüfstand  gemessenen  Anregungs¬ 
kräften  eines  Hinterachsgetriebes  am  Prüfstand  (mitte) 
und  der  gemessenen  FRF  des  Fahrzeugs  (oben) 


Dieses  Ergebnis  zeigt,  dass  die  Kraft  f  aus  der  Mes¬ 
sung  der  Koppelkraft  fs  sowie  der  zugehörigen  Ver¬ 
schiebung  u2  des  Lagerpunktes  am  (nachgiebigen) 
Prüfstand  bestimmt  werden  kann.  Für  die  Berechnung 
wird  darüberhinaus  nur  die  dynamische  Nachgiebig¬ 
keit  des  Getriebes  am  Anbindungspunkt  benötigt. 
Diese  kann  sowohl  mittels  Shaker  oder  Impulsham¬ 
mermethode  gemessen,  als  auch  mittels  FE  Berech¬ 
nung  ermittelt  werden.  Sowohl  die  Eigenschaften  des 
Prüfstandes  als  auch  die  anregenden  Zahnkräfte  müs¬ 
sen  hingegen  nicht  bekannt  sein. 
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